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This  report  was  prepared  by  Simula  Inc.  under  Contract  No.  DTFA03-83-C-00036 
with  the  Federal  Aviation  Administration  (FAA)  Technical  Center,  where  L.  M. 
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development  and  validation  of  program  SOM-TA,  and  the  preparation  of  this  re¬ 
port.  Data  for  model  validation  were  provided  by  the  Protection  and  Survival 
Laboratory  of  the  FAA  Civil  Aeromedical  Institute,  where  the  tests  were  con¬ 
ducted  under  the  leadership  of  Mr.  R.  F.  Chandler. 
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EXECUTIVE  SUMMARY 


Program  SOI i- T A  (Seat/Occupant  Model  -  Transport  Aircraft)  has  been  developed 
for  use  in  evaluating  the  crashworthiness  of  transport  aircraft  seats  and  re¬ 
straint  systems.  It  combines  a  three-dimensional  dynamic  model  of  the  human 
body  with  ?  finite  element  model  of  the  seat  structure.  The  program  allows 
simulation  of  one,  two,  or  three  passengers  of  the  same  or  different  sizes. 

It  is  intended  to  provide  the  design  engineer  a  tool  with  which  he  can  analyze 
the  structural  elements  of  the  seat  as  well  as  evaluate  the  dynamic  response 
of  the  passengers  during  a  crash. 

The  program  is  documented  in  two  volumes.  Volume  I  discusses  development  of 
the  seat  and  occupant  models,  validation,  and  organization  of  the  computer 
program.  Volume  II  contains  instructions  for  preparing  input  data  and  operating 
the  program. 

The  occupant  model  consists  of  12  masses  that  represent  the  upper  arid  lower- 
torso,  neck,  head,  and  two  segments  for  each  of  the  arms  and  legs.  An  optional 
model  of  the  human  body  includes  beam  elements  in  the  spine  and  neck,  but  is 
restricted  to  two-dimensional  motion. 

External  forces  are  applied  to  the  occupant  by  the  cushions,  floor,  and  re¬ 
straint  system.  Interface  between  the  occupant  arid  seat  is  provided  by  the 
seat  bottom  cushion,  back  cushion,  and  an  optional  headrest.  Contact  between 
the  occupant(s)  and  the  back  of  the  seat  row  in  front  can  also  be  simulated. 

The  restraint  system  can  consist  of  a  lap  belt  alone  or  combined  with  a  single 
shoulder  belt,  ovev  either  shoulder,  or  a  double-strap  shoulder  harness.  Each 
component  of  the  restraint  system  can  be  attached  to  either  the  seat  or  the 
aircraft  structure. 

The  seat  structure  is  modeled  using  the  finite  element  method  of  analysis, 
selected  because  it  is  not  dependent  on  previous  testing,  and  it  has  the  flex¬ 
ibility  to  deal  with  a  wide  range  of  design  concepts.  The  SOM-TA  seat  analysis 
uses  three-dimensional  beam  elements.  It  has  the  capability  to  model  large 
displacements,  nonlinear  material  behavior,  local  buckling,  and  various  end 
release  conditions  for  beam  elements. 

The  digital  computer  program  based  on  the  occupant  and  seat  models  described 
above  has  boon  written  entirely  in  FURTRAH  to  ensure  a  high  degree  of  compati¬ 
bility  with  various  computer  systems.  During  development,  the  program  has  been 
run  on  Ilil-i  and  Digital  computer  systems.  Output  data  include  occupant  segment 
positions,  velocities,  and  accelerations;  restraint  system  and  cushion  forces; 
injury  criteria;  anu  details  of  contact  between  the  occupant  and  the  seat  in 
front.  Seat  output  includes  nodal  displacements,  element  stresses,  and  forces 
at  the  points  of  attachment  t.u  the  aircraft  structure. 

Dynamic  tests  of  transport  seats  conducted  by  the  Federal  Aviation  Administra¬ 
tion  livil  Aerouedical  institute  in  1983  were  used  in  validation  of  Program 
Son-'IA.  Three  grouuction  seat  designs  and  various  modifications  to  the  produc¬ 
tion  design  were  tested.  The  design  that  was  chosen  for  validation  utilized  a 
production  seat  modi  tied  by  incorporating  energy-absorption  capability  in  the 
longitudinal  direction.  This  seat  produced  more  interesting  results  for  vali¬ 
dation  than  others  because  it  underwent  large  deformations  while  maintaining 
structural  integrity  in  the  dynamic  test. 


1.0  INTRODUCTION 


The  design  of  crashworthy  seats  and  restraint  systems  for  transport  aircraft 
presents  a  complex  engineering  problem,  the  solution  of  which  can  be  greatly 
aided  by  sufficiently  rigorous  analytical  techniques.  The  crash  environment 
can  vary  widely  from  one  accident  to  another,  so  a  great  number  of  condi¬ 
tions  must  be  evaluated  to  establish  those  critical  to  occupant  survival.  A 
crashworthy  seat  should  include  the  capacity  to  absorb  energy  through  con¬ 
trolled  deformation,  thus  reducing  the  accompanying  loads  (reference  1). 

In  the  initial  design  phases,  it  is  desirable  to  evaluate,  in  some  detail, 
existing  seats  and  restraint  systems  in  their  surroundings,  thus  establishing 
existing  weaknesses,  it  is  then  desirable  to  make  modifications  and  to  evalu¬ 
ate  the  effect  of  these  modifications  on  improving  the  performance  of  the 
system.  These  evaluations  must  be  conducted  for  a  yreat  many  of  the  possible 
crash  environrents,  thus  constituting  a  relatively  large  matrix.  Testing  is 
costly  and  is  often  time  consuming  since  design  modifications  must  be  developed 
and  fabricated  prior  to  testing.  Therefore,  an  analytical  technique,  such  as 
this  effort,  is  cost  effective. 

A  number  of  one-,  two-,  and  three-dimensional  mathematical  models  of  the  human 
body  have  been  developed  for  crash  survivability  analysis.  These  models  vary 
in  complexity  and  possess  from  one  to  forty  deyrees-of-freedon.  The  simplest 
models  have  been  developed  primarily  for  prediction  of  injury  to  a  single 
component  or  subsystem  of  the  body,  such  as  vertebral  injury  due  to  the  verti¬ 
cal  force  that  might  be  experienced  in  a  crash  or  the  firing  of  an  ejection 
seat.  Simulation  of  the  three-dimensional  response  of  the  entire  body  re¬ 
quires  many  more  degrees-cf-freedoiu,  but  permits  more  general  use.  host  of 
the  three-dimensional  models  have  been  developed  for  use  in  evaluation  of 
automobile  interior  designs  with  respect  to  injuries  caused  by  secondary  im¬ 
pacts,  such  as  the  three-dimensional  models  described  in  references  1  through 
b.  Seats  have  been  represented  in  a  very  simple  manner  because  in  automobiles 
the  role  of  the  seat  design  in  determining  occupant  survival  is  minimal. 
Therefore,  a  simulation  model  intended  specifically  for  the  aircraft  applica¬ 
tion  is  desired. 
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The  development  of  a  three-dimensional  mathematical  model  of  a  transport  air¬ 
craft  seat,  including  passengers,  and  restraint  system  is  described  in  this 
report.  This  model  forms  the  basis  for  a  simulation  computer  program  that  has 
been  written  specifically  for  use  in  crashworthy  design  and  analysis  of  trans¬ 
port  aircraft  seats  and  restraint  systems,  Program  SOM-TA.  The  program,  which 
combines  a  finite  element  model  of  the  seat  structure  with  lumped -parameter 
models  of  the  aircraft  passengers,  has  been  organized  so  as  to  minimize  the 
volume  and  complexity  of  input  data  and  to  focus  on  seat  and  restraint  system 
design  parameters. 

Proyraiii  tiUM-TA  seat  and  occupant  models  are  based  on  those  presently  used  in 
Program  bOh- LA  (Seat/Occupant  Model  -  Liyht  Aircraft),  (reference  6).  The 
proyraia  was  modified  to  allow  simulation  of  one,  two,  or  three  passengers  of 
the  same  or  different  sizes.  The  capacity  of  the  seat  model  has  been  increased 
to  accommodate  the  more  complex  transport  seat  structures.  The  seat  model  has 
also  been  modified  to  simulate  floor-warping  conditions. 

In  this  report,  the  details  of  the  occupant  and  seat  models  are  presented, 
alony  with  the  results  of  model  validation.  Operation  of  the  computer  pro¬ 
gram  is  described  in  Volume  II,  Program  SOM-TA  User  Manual. 
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2.0  OCCUPANT  MODEL 


The  SOM-TA  mathematical  model  includes  a  lumped-parameter  representation  of 
the  aircraft  occupant(s)  and  a  finite  element  characterization  of  the  seat. 
Interface  between  the  seat  and  its  occupant  is  provided  by  seat  cushions  and  a 
restraint  system,  which  consists  of  a  lap  belt  and,  if  desired,  a  single-strap 
or  double-strap  shoulder  harness.  The  response  of  the  occpant(s)  and  seat 
can  be  predicted  for  any  given  set  of  aircraft  impact  conditions,  including 
the  initial  velocity,  attitude  and  the  input  acceleration. 

As  mentioned  previously,  both  a  three-dimensional  occupant  model  consisting  of 
rigid  links  and  a  two-dimensional  model  with  deformable  spinal  elements  are 
included.  This  chapter  provides  a  discussion  of  the  development  of  both  oc¬ 
cupant  models,  including  details  of  the  approach  used  to  formulate  the  equa¬ 
tions  of  motion  and  the  techniques  used  for  their  solution.  Section  2.1  pre¬ 
sents  the  three-dimensional  model  and  section  2.2,  the  two-dimensional  model. 
Subsequent  sections  cover  aspects  of  the  equations  which  apply  to  both  models, 
specifically  the  body  joint  model,  the  treatment  of  external  forces,  and  body 
dimensions  and  properties. 

2.1  THREE-DIMENSIONAL  OCCUPANT  MODEL 


The  three-dimensional  mathematical  model  of  the  aircraft  occupant  is  made  up 
of  twelve  rigid  segments,  as  shown  in  figure  1.  This  number  is  thought  to  re 
present  the  minimum  segments  that  will  permit  accurate,  meaningful  simulation 
a  of  three-dimensional  response.  A  greater  number  of  segments  might  possibly 
improve  the  accuracy  of  simulation  but  would,  in  turn,  also  increase  program 
execution  cost.  Arm  and  leg  segments  are  included  to  enable  prediction  of  in 
juries  to  these  extremities.  Although  leg  and  arm  injuries,  in  themselves, 
may  not  be  as  serious  as  head  or  chest  injuries,  they  may  prevent  egress  from 
a  stricken  aircraft  and  increase  the  potential  hazard  to  the  occupants  of  a 
postcrash  fire. 

Each  of  the  body  joints,  with  the  exception  of  the  elbow,  knee,  arid  head-neck 
joints,  possesses  three  rotational  degrees  of  freedom.  Because  of  the  hinge- 
type  motion  of  a  forearm  or  lower  leg  relative  to  an  upper  arm  or  thigh,  the 
position  of  each  of  these  segments  is  described  by  one  additional  angular 


■WVCi&g Si 


» » *»>  aVvv  VkLv 


/ 


©SEGMENT  MASS  CENTER 
O  JOINT  ®,  @ . © 


Figure  1.  Twelve-segment  ( three-dimensional ) 
occupant  model. 


coordinate,  as  also  is  the  position  of  the  head  relative  to  the  neck.  There¬ 
fore,  the  occupant  system  possesses  a  total  of  29  degree s-of- freedom. 

2.1.1  Coordinate  Systems 

Fixed  at  the  center  of  mass  of  each  of  the  12  segments  is  a  right-handed  Car¬ 
tesian  coordinate  system.  For  segment  n  (n  =  1,  2,...,  12)  the  local  coordi¬ 
nate  system  is  denoted  by  axes  (x  ,  v  ,  z  ).  Positive  directions  are  defined 
•7  n  n  n 

such  that  when  the  body  is  seated  as  shown  in  figure  2,  with  the  torso  and  head 


upright,  the  upper  arms  parallel  to  the  torso,  and  the  elbows  and  knees  bent 
at  right  angles,  positive  xn  is  directed  forward,  yp  to  the  left,  and  zn  upward. 


Maty 

In  order  to  describe  a  general  position  of  the  body,  it  is  necessary  to  relate 
the  orientation  of  each  segment  (x  ,  yn ,  z  )  to  the  Inertial  system  (X,  Y,  Z). 

The  angular  relationship  between  the  local,  segment-fixed  coordinates  and  the 
inertial  system  can  be  expressed  by  the  transformation 


(1) 


Because  three  angular  coordinates  can  be  used  to  define  the  rotation  of  a  given 
segment,  it  is  convenient  to  utilize  a  set  of  coordinates  that  will  suffice  as 
generalized  coordinates  in  the  formulation  of  the  equations  of  motion.  A  system 
of  Eulerian  angles  provides  a  convenient  set  of  three  independent  angular  co¬ 
ordinates.  Assuming  that  the  local  (xn,  yn,  zn)  system  is  initially  coinci¬ 
dent  with  the  inertial  (X,  Y,  Z)  system,  the  Euler  angles  are  a  series  of  three 
rotations,  which,  when  performed  in  the  proper  sequence,  permit  the  system  to 
attain  any  orientation  and  uniquely  define  that  position.  The  particular  set 
of  Euler  angles  selected  for  use  here  is  illustrated  in  figure  3  and  defined 
as  follows: 


1.  A  positive  rotation  ip  about  the  Z-axis,  resulting  in  the  primed  (x ' , 
y ' ,  z 1 )  system. 

2.  A  positive  rotation  0  about  the  y'-axis  resulting  in  the  double-primed 
(x",  y",  z")  system. 

3.  A  positive  rotation  about  the  x"-axis  resulting  in  the  final  (x, 
y,  z)  system. 


■KU/! 

Hv 


6 


.•V-.V 


Figure  3.  The  Euler  angles. 


In  order  to  determine  the  elements  of  the  transformation  matrix  [Tn],  it  is 
necessary  to  consider  the  matrix  equations  that  indicate  the  three  individual 
rotations  previously  described  above.  Referring  again  to  these  definitions 
of  \(j,  e,  and  <j>,  the  following  equations  are  obtained: 


cos\p  -sim})  0 
si n\p  cos\Jj  0 


0  1 


cose  0  sine 
0  1  0 


-sine  0  cose 


1  0  0 
0  cos 4>  -sin4> 
0  sine))  cos 4) 


Writing  equations  (2)  through  (4)  in  abbreviated  form 


{R}  =  O]  {r*  } 


(r ' }  =  [e]  (>""  > 


(r"}  =  [$]  {r} 


{R}  -  !>]  [6]  [<t>]  {r}  (5) 

where  {R}  represents  the  components  of  a  vector  in  the  inertial  system  and  {r} 
represents  the  same  vector  in  the  final  (x,  y,  z)  system.  Performing  the  matrix 
multiplications  Indicated  in  equation  (5),  the  elements  of  the  transformation 
matrix  in  equation  (1)  are  obtained: 
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The  constraint  of  hinge-type  rotation,  at  the  elbows  and  knees,  requires  the 
use  of  one  additional  angular  coordinate  to  define  the  position  of  each  of  the 
forearm  and  lower  leg  segments.  Referring  to  figure  4,  the  angular  position 
of  the  forearm  segments  (  n  =  5,  7)  is  given  by 


r  \ 


sin 

0 

cos 

aH 

1 

f  -\ 

Tn 

0 

1 

0 

I 

-cos  a ^ 

0 

si  n 

aH 

ZJ 

-  ^ 

M 

L  J 

(7) 


and  the  lower  leg  segments  (m  =  9,  11)  by 
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a  cos  cos  0n  sinci^  -  cos  ^  sin0n  cos  4>  cosa£ 

-  sin  i>  sin  <t)  cos 

l 

T12  b  cos  iiin  sin  9n  sin  4>  -  sin  \|i  cos  <J» 

T13  s  cos  4>n  cos  9n  cosa^  +  cos  \|>n  sin  0n  cos  <j>n  sina  £ 
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a 
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T32  -cosens1M,n 
£ 

T33  =  -  sin  0n  cosa^  +  C0o3  cost})  n  sinoi^  (9) 


From  equations  (6)  and  (8),  the  elements  of  the  transformation  matrix  elements 
For  the  legs  (segments  9  and  11)  are  obtained: 
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(10) 


Segment  12,  the  neck,  was  added  later.  The  hinge-type  joint  between  the  neck 
and  the  head  is  treated  similarly  to  the  knee  joints,  adding  the  generalized 
coordinate  012,  as  illustrated  in  figure  4. 


Having  developed  the  relationships  expressed  in  the  equations  (1)  through  (10), 
the  position  of  the  occupant  can  be  described  by  the  following  set  of  general¬ 
ized  coordinates: 
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(11) 

The  above  coordinates  include  the  Cartesian  coordinates  of  the  mass  center  of 
segment  1  (Xj,  Y^,  Z^),  selected  as  a  reference  point  on  the  body,  seven  sets 
of  Eulerian  angles,  and  the  five  additional  angular  coordinates  for  the  elbows, 
knees,  and  neck.  Positions  of  the  segment  mass  centers  are  presented  in  appen¬ 
dix  A. 


2.1.2  Lagrange's  Equations 

The  response  of  the  occupant  system  is  described  by  Lagrange's  equations  of 
motion,  which  are  written  for  the  29  generalized  coordinates.  The  equations 
are  developed  according  to 


iL  (iL)  _  =  n . 

dt  3q0  *J  J 


(j  =  1,  2,...,  29) 


where  Lis  the  Lagrangian  function 

L  ■  T  -  V  (13) 

t  represents  time,  Qj  are  the  generalized  forces  not  derivable  from  a  potential 
function.  (Forces  that  are  derivable  from  a  potential  function  are  obtained 
from  L,  and  T  and  V  are  the  system  kinetic  and  potential  energies,  respectively.) 


Because  the  system  being  treated  does  not  involve  any  velocity-dependent  poten¬ 
tials,  equation  (12)  can  be  written  as 


d  ,3T  s  3T  ,  31 
iT  117“  > 


U  /  U  I  \  V  » 

dt  "  3qj  3qj 


Q,i 


(j  -  1,  £ .  29) 


(14) 


The  system  kinetic  energy  contains  both  translational  and  rotational  parts: 
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t4  i 
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A  y  (i  to  2  +  i  to  2  +  i  u  2) 

2  n-1  xn  xn  yn  zn  zn 


(15) 


where  M  is  the  mass  of  segment  n  and  I  ,  I  ,  and  I  are  toss  moments  of 
n  xn  yn  n 

inertia  of  segment  n  with  respect  to  the  local  coordinate  axes  (xn,  yn>  zn), 

assumed  to  be  principal  moments  of  inertia. 


The  absolute  velocities  of  the  12  mass  segments  required  for  the  translational 
kinetic  eneryy  must,  of  course,  be  written  as  runctions  of  the  generalized  co¬ 
ordinates  and  generalized  velocities  in  order  to  use  equation  (14).  The  angu¬ 


lar  velocity  components  (to  ,  <o  ,  to  )  seen  in  equation  (15)  are  parallel  to 

n  yn  n 


the  local  (x  ,  yp ,  zn)  coordinate  systems.  These  angular  velocity  components 


cannot  be  used  directly  in  Lagrange's  equations  because  they  do  not  correspond 
to  the  time  derivatives  of  any  set  of  coordinates  that  specify  the  position  of 
the  segment.  They  must  be  written  as  functions  of  the  generalized  coordinates, 
using  the  generalized  angular  velocities  (j>  ,  0  ,  dp)»  which  are  parallel  to 


the  axes  Z,  y^,  and  xj|,  respectively. 


An  arbitrary  angular  velocity  of  segment  n,  w  ,  can  be  expressed  as  a  function 
of  the  generalized  angular  velocities  according  to 


O)  =  ^  +0  f  <t> 

-n  -n  -n  -n 


(16) 


Referring  to  figure  3,  0n>  and  j>n  do  not,  in  general,  form  a  mutually  per¬ 

pendicular  vector  triad,  (j;  and  j>  are  both  perpendicular  to  <3  but  are  not  nec¬ 
essarily  perpendicular  to  each  other.)  However,  they  can  be  considered  as  a 
nonorthogonal  set  components  of  w  since  their  vector  sum  is  equal  to  w. 
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Summing  the  orthogonal  projections  of  Jn,  ®n»  and  tn  on  the  ^xn’  yn*  zn^  dXes 
yields  the  angular  velocity  components  required  for  the  kinetic  energy  expres¬ 


sion: 


“x  ■  K  ■  K  5<"  K 
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“yn  ■  *„  cos  9n  s1n  *n  +  9n  cos  *n 


“j  "  "ii  cos  6n  cos*n  min  s1n*n 
n 


The  system  potential  energy  Is  simply  gravitational  potential,  which  is  writ¬ 


ten  as 


v  -  f  (\,g  (Z„  -  z  )  (is; 

n=l  o 

^jgjjp  where  g  is  the  acceleration  due  to  gravity  and  Zn  is  an  arbitrary  datum. 


2.1.3  Matrix  Equations 

For  the  purposes  of  computation,  the  equations  of  motion  are  rewritten  in  the 
following  fo rtn: 

[A(q)]  {q>  =  (B(q,q) }  +  {P(q)>  +  {R(q,q) >  +  (Q(q.q) >  (19) 

where  the  elements  of  the  inertia  matrix  [A]  and  the  vector  {B}  are  derived 
from  the  kinetic  energy  derivatives  of  Lagrange's  equations.  In  other  words. 


Vvs 


d _  /  3T  .  aT 

dt  Lqj  '  3<>j 


-  Bj  ( q  i »  ^2 » •  •  * »  9  29 »  c^2,••*,  ^29^ 


(j  =  1,  2 . 29) 


The  force  vector  { P }  is  derived  from  the  system  potential  energy  according  to 
Pj  (qr  q2,....  q29)  =  -  ||j  (j  =  1,  2,....  29)  (21) 

Both  {K}  and  {q}  are  vectors  of  generalized  forces  derived  from  the  right-hand 
side  of  Lagrange's  equations.  The  vector  { R }  describes  the  resistance  of  the 
body  joints  to  rotation,  discussed  in  detail  in  section  2.3;  { Q }  is  the  vector 
of  generalized  external  forces,  discussed  in  detail  in  section  2.4. 

2*2  TWO-DIMENSIONAL  OCCUPANT  MODEL 

In  order  to  achieve  economical  program  solutions  for  cases  where  occupant  re¬ 
sponse  is  expected  to  be  symmetrical  with  respect  to  the  X-Z  plane,  a  two- 
dimensional  occupant  model  option  was  included  in  Program  SOM-TA.  This  two- 
dimensional  model  was  configured  like  the  three-dimensional  model  shown  in 
figure  1,  with  the  exception  of  all  joints  being  hinge-type  joints.  Because 
of  the  potential  for  vertebral  injury  in  aircraft  accidents  that  involve  a 
significant  vertical  component  of  impact  velocity,  some  measure  of  vertebral 
loading  was  considered  desirable  in  the  occupant  model.  The  two-dimensional 
occupant  model  was  configured  to  include  beam  elements  in  both  the  torso  and 
neck,  as  shown  in  figure  5,  replacing  joints  that  exist  in  the  three-dimensional 
model.  The  two-dimensional  model  has  11  degrees-of- freedom,  as  illustrated  in 
figure  6.  Simulation  of  purely  vertical  (+G  )  impact  with  the  three-dimensional 
occupant  model  produced  less  bending  of  the  torso  than  observed  in  tests  with 
either  dummies  or  human  cadavers.  Therefore,  in  order  to  produce  the  moments 
on  the  vertebral  column  that  are  induced  by  +G  acceleration  in  both  human  and 
dummy  occupants,  the  mass  centers  of  the  torso  segments,  segments  1  and  2,  are 
offset  fran  the  vertebral  column  by  the  dimensions  e^  and  e,,,  respectively. 

The  dimension  e^,  by  which  the  center  of  mass  of  the  head  segment  is  placed 
forward  of  the  neck,  is  used  in  both  three-dimensional  and  two-dimensional 
model s. 

Development  of  the  equations  of  motion  for  the  two-dimensional  occupant  follows 
a  procedure  similar  to  that  described  in  section  2.1  for  the  three-dimensional 
model.  However,  the  procedure  is  simplified  because  the  transformation  between 


A  detailed  representation  of  the  kinematics  and  mechanics  of  the  spinal  column 
would  be  complicated.  The  presence  of  the  Intervertebral  disks  leads  to  high 
flexibility  in  bending  and  high  stiffness  in  compression.  In  addition,  the 
column  is  tapered  and  possesses  considerable  initial  curvature.  A  simplified 
model  is  incorporated  into  SQM-TA  by  using  continuous  beam  elements  for  the 
neck  and  spine.  These  beam  elements  are  intended  to  model  the  flexural,  as 
well  as  axial,  motion  of  the  spine  and  neck  and  are  subject  to  the  following 
assumptions.  The  deformed  beam  elements  take  the  shape  of  circular  arc  seg¬ 
ments,  therefore  assuming  flexure  to  be  primarily  due  to  the  applied  bending 
moment.  This  Is  equivalent  to  the  assumption  in  beam  theory  that  the  span/ 
depth  ratio  is  large  and  that  the  modulus  of  elasticity  and  shear  modulus  are 
of  the  same  order  of  magnitude  so  that  shear  deformation  is  negligible. 


As  illustrated  in  figure  7,  The  position  of  one  end  relative  to  the  other  In 
terms  of  arc  length  and  angular  coordinates  at  each  end  of  the  arc,  is  then 
given  by 


Positions  of  the  mass  centers  for  the  eleven  body  segments  of  the  two-dimensional 
model  are  presented  in  appendix  B. 

2.3  JOINT  RESISTANCE 

The  form  of  the  joint  resistance  vector  { R}  in  equation  (19)  depends  on  the 
user's  selection  of  occupant  type  -  either  dummy  or  human.  Although  both  joint 
models  contain  the  same  types  of  elements,  a  nonlinear  torsional  spring  and  a 
viscous  torsional  damper,  the  relative  contributions  of  each  of  these  elements 
determine  the  type  of  occupant. 


The  11  body  joints  for  the  three-dimensional  model,  illustrated  in  figure  1, 
are  defined  as  fol lows: 
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Figure  7.  Position  of  segment  Z  relative 
to  segment  1. 


Joint  1  -  Back,  between  12th  thoracic  and  1st  lumbar  vertebrae 

Joint  2  -  Torso-neck,  between  7th  cervical  and  1st  thoracic  vertebrae 

Joint  3  -  Right  shoulder 

Joint  4  -  Right  el  bow 

Joint  5  -  Left  shoulder 

Joint  6  -  Left  elbow 

Joint  7  -  Right  hi p 

Joint  8  -  Right  knee 


Joint  9  -  Left  hip 
Joint  10  -  Left  knee 

Joint  11  -  Head-neck,  at  occipital  condyles 

The  angular  displacement  of  joint  1  from  its  reference  position  (figure  2)  is 
given  by  B,.  If  (i„.  Jn.  *m)  and  U„_  d„.  kn)  are  triads  of  unit  vectors  In 
the  local  coordinate  systems  of  two  adjacent  segments  connected  at  joint  1,  as 
shown  in  figure  8,  the  joint  angle  Is  given  by 

Si  -  cos-1  {km  •  kn)  (23) 

where  (k  •  fe  )  is  the  scalar  product.  Considering  the  geometry  of  the  occu¬ 
pant  model  in  the  reference  position,  the  8^  for  the  11  joints  are  given  by 
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Figure  8.  Joint  angle  3^  between  segments  m  and  n. 


If  at  each  joint  i,  a  moment  and  a  torsional  clamper  with  coefficient  act 


to  resist  motion  of  the  joint,  then  the  virtual  work  done  on  the  systan  as  each 
joint  i  undergoes  a  virtual  displacement  68^  is 
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6W  =  -  l  (N1  68,  +  d1  B1  681 ) 


(25) 


Since  the  3.  are  functions  of  the  generalized  coordinates  ,  the  virtual  dis¬ 
placements  68.  can  be  expressed  in  terms  of  corresponding  virtual  displacements 


of  the  q..  In  general,  such  an  expression  would  take  the  form 
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where  the  partial  derivatives  3^/aqj  are  functions  of  the  generalized  coordi¬ 
nates.  Substituting  into  equation  (25)  gives 


6W  *  -  \  f  (M,  +  Ji  V 

1  =  1  J  ’"‘1  J 


Changing  the  order  of  summation,  equation  (27)  can  be  written  in  the  general 


6  W  =  l  Rj  5qj  (28) 

j=l  J  J 

where  R.  are  the  generalized  forces  acting  on  the  system.  As  seen  in  equation 
(19)»  the  generalized  forces  are  treated  as  two  distinct  types t  joint  resist¬ 
ance  forces  and  external  forces.  Since  the  joint  resistance  terms  are  being 
treated  here,  the  generalized  joint  forces  referred  to  as  Rj  will  be  consid¬ 
ered  alone.  Equation  (28)  becomes,  more  specifically 


rj  ^ 


and,  from  equation  (27),  Rj  can  be  written 


■- 


11  . 


(j-1,  2,...,  29) 


As  mentioned  earlier  in  this  section,  the  type  of  occupant  is  determined  by 
the  relative  contributions  of  and  J1  to  the  Rj  terms.  For  the  dumn\y  joint, 
the  resisting  torque  M.  is  constant  throughout  the  normal  range  of  joint  mo¬ 
tion  and  increases  rapidly  along  a  third-order  curve  to  a  higher  value  at  the 
limiting  displacement  Bs  ,  as  shown  in  figure  9.  The  normal  values  MD  are 

set  equal  to  those  resulting  from  the  joint-tightening  procedure  of  SAE  Recom¬ 
mended  Practice,  Anthropomorphic  Test  Device  for  Dynamic  Testing  -  SAE  J963. 
That  is,  the  body  joints  will  just  support  a  1-6  load  in  the  reference  (seated) 


u.  uv.v  -.rw\  ".ur.iv.  i  vr.v.y-.i 


Angular  displacement  -  0^ 


Figure  9.  Dumny  joint  resisting  torque. 


position,  with  the  exception  of  the  torso  joints,  which  will  support  a  2-G  load. 
In  addition  to  M^,  a  small  viscous  damping  term  with  constant  Is  Included 
for  energy  dissipation. 

The  resistance  of  each  human  joint  consists  of  up  to  three  terms.  The  primary 
resisting  force  during  normal  joint  rotation  is  a  viscous  damping  term  with 
constant  coefficient  J..  In  a  manner  similar  to  the  case  of  the  dummy,  a  re¬ 
sisting  torque  is  applied  at  the  limit  of  the  joint  range  of  motion,  as  shown 
in  figure  10(a).  An  additional  term  used  to  simulate  muscle  tone  is  the  mo¬ 
ment  M 1 ,  which  drops  to  zero  after  a  small  angular  displacement  from  the  ini¬ 
tial  position,  provided  that  the  crash  deceleration  is  sufficient  to  overcome 
1 t  (figure  10(b) ) . 

For  the  two-dimensional  model,  the  equations  presented  in  this  section  are 
correct  except  that  summations  are  performed  for  8  joints  and  11  degrees  of 
freedom,  rather  than  the  11  and  29  presented  here,  respectively. 


Figure  10.  Human  joint  resisting  torques:  (a)  displacement 
limiting  moment  and  (b)  muscular  resistance. 


2.4  EXTERNAL  FORCES 

The  vector  of  gereralized  external  forces  {Q}  is  developed  in  a  manner  similar 
to  that  discusseo  in  the  previous  section  for  the  joint  resistance  vector. 
Equations  presented  in  this  section  include  summations  whose  upper  limits  are 
correct  for  the  three-dimensional  model,  with  12  segments  and  29  degree  s-of- 
freedotn.  The  approach  for  the  two-dimensional  model  is  the  same,  except  for 
the  use  of  11  segments  and  11  degrees- of- freedom. 

The  resultant  external  force  acting  on  segment  i  is  given  by 

F,  =  F  i  +  Fv  0  +  F  k  (31) 

-i  x.j  *  i  L\ 


where  Fy  ,  Fv  ,  and  F7  are  components  in  the  inertial  (X,  Y,  Z)  system.  The 
Ai  ri  4i 

absolute  position  of  the  point  P ^  on  'gment  i,  where  the  resultant  force  acts, 
can  be  represented  by 
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(32} 


,1 


VV+V+V  (32) 

As  the  resultant  force  applied  to  each  segment  1  undergoes  a  virtual  displace¬ 
ment  <srp^,  having  components  (6Xp^,  SYp  ,  cZp  ),  the  virtual  work  on  the  system 

done  by  the  is 


5W  =  I  (F  SXp  +  Fy  6Yp  +  F,  6Zp  ) 
i=l  xi  *1  '1  pi  L\  pi 


Writing  the  virtual  displacement  components  in  terms  of  the  generalized  coor¬ 
dinates  q.: 

J 


29  3XP 

6XP.  *  J.  ‘SqT’  \ 
i  jal  J 


29  3^P 


29  32P, 


i  j=l 


aqr0Mj 


results  in 


29  12 


>i  iJ=i 


Using  equation  (28) 


-  l  Q.  6q. 
j=l  J  J 


yields  the  components  of  the  generalized  external  force  vector: 
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(F 
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3Xp  3Yp 

_ i  i  x  P 

W  % 


(3G) 


The  external  forces  acting  on  the  body  segments  can  be  characterized  as  either 
contact  forces  or  restraint  forces.  These  forces  are  discussed  in  further  de¬ 
tail  in  the  sections  following. 

2.4.1  Contact  Forces 


The  contact  forces  applied  to  the  occupant  in  all  cases  are  those  forces  exerted 
by  the  cushions  and  floor,  illustrated  in  figure  11.  The  forces  of  the  seat 
back  pass  through  the  mass  centers  of  the  upper  torso  and  head  segments,  and  the 
normal  forces  applied  to  the  legs  by  the  seat  bottom  cushion  pass  through  the 
mass  centers  of  the  thigh  segments.  The  contact  surface  for  the  lower  torso  is 
not  located  at  the  mass  center  of  that  segment  but  is  an  ellipsoid  whose  major 
axis  passes  through  both  hip  joints.  The  force  of  the  back  cushion  and  the  normal 
component  of  the  seat  bottom  cushion  force  both  pass  through  the  center  of  this 
ellipsoid,  i.e.,  through  a  point  midway  between  the  hips.  The  normal  components 
of  the  floor  and  footrest  forces  are  applied  to  the  body  at  the  lower  ends  of 
the  leg  segments.  If  contact  between  the  occupants  and  the  seat  in  front  of 
them  is  to  be  simulated,  eight  surfaces  on  each  seat  back  are  defined  by  user 
input.  If  contact  occurs  between  an  occupant  and  a  seat  back,  forces  are 
computed  in  the  same  manner  as  for  cushions. 

Fach  of  the  forces  mentioned  above  acts  normal  to  the  surface  applying  the 
force,  or,  in  other  words,  its  direction  is  determined  by  the  plane  of  the  sur¬ 
face,  whether  cushion  or  floor.  As  shown  in  figure  11,  friction  forces  are 
also  applied  by  the  seat  bottom  cushion  and  the  floor.  The  friction  force  is 
computed  as  the  product  of  the  coefficient  multiplied  by  the  normal  force.  It 
is  applied  in  a  direction  opposite  to  the  tangential  component  of  relative  ve¬ 
locity  between  the  occupant  segment  and  the  appropriate  cushion  or  floor  sur¬ 
face.  In  order  to  avoid  abrupt  changes  in  direction  of  the  friction  force, 
the  force  is  reduced  sinusoidally  when  the  magnitude  of  the  tangential  veloc¬ 
ity  drops  below  a  predetermined  limit.  (A  limiting  value  of  1  ft/sec  has  pro¬ 
duced  satisfactory  results  in  SOM-TA  and  is  thus  used  by  the  program.) 
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Figure  11.  External  forces  of  cushions  and  floor. 

All  contact  forces  are  calculated  by  first  determining  the  penetration  of  a 
contact  surface  on  the  occupant  into  a  surface  with  known  force-deflection 
characteri sties.  lln'rig  the  seat  cushion  force  as  an  example,  the  pertinent 
dimensions  of  the  seat  and  the  parameters  required  to  determine  the  penetra¬ 
tion  of  the  abdanino-pel vie  segment  (segment  1)  into  the  cushion  are  illus¬ 
trated  in  figure  12.  Xp  and  Zp  are  coordinates  of  the  cei.ter  of  the  contact 
surface  of  segment  1,  and  R^  is  the  radius  of  the  contact  surface  in  the  ( 
z^)  plane.  (Although  this  contact  surface  is  an  ellipsoid,  cross-sections 
parallel  to  the  (x^  -  z^)  plane  are  circular.  The  dimensions  of  the  contact 


Figure  12.  Seat  cushion  deflection. 

surfaces  will  be  discussed  in  section  2.5.)  The  position  of  the  seat  pan  is 

defined  by  its  height  Z <.  above  the  origin  of  the  aircraft  coordinate  system  and 

the  angle  e,.  that  it  makes  with  the  aircraft  (X  -  Y  )  plane.  The  unloaded 

thickness  of  the  seat  cushion  is  t  ,  and  the  loaded  thickness  under  segment  1 

6 

is  t.  Summing  the  dimensions  in  the  Z A  direction  gives 

Zp  ~  Z(.  +  (R1  +  t)/cos  6S  +  Xp  tan  @s  (37) 

Solving  equation  (37)  for  the  cushion  thickness, 

t  =  (Z..  -  Z  )  cose  -  R,  -  X-  sine  (38) 

F  S  S  1  P  S 

The  deflection  of  the  seat  cushion  is  then 

6  =  t  -  t  (39) 

c  e 

and  the  force,  which  is  assumed  to  act  normal  to  the  plane  of  the  seat  pan  and 
pass  through  the  center  of  curvature  of  the  contact  surface,  is  calculated  from 
deflection  according  to 


82  C I  0  0  3 


To  each  normal  force,  a  damping  term  Is  applied  which  is  proportional  to  the 
deflection  rate.  The  damping  coefficient  is  based  on  the  Rayleigh  formulation 
in  which  the  coefficient  is  proportional  to  both  mass  and  stiffness  according 
to 


C  ®  2  Bni  +  a K  (41) 

For  a  multidegree-of-freedom  system  there  will  be  a  discrete  damping  coeffi¬ 
cient  associated  with  each  characteristic  mode.  In  a  continuous  system  there 
will  be  an  infinite  number  of  coefficients,  although  several  modes  will  gener¬ 
ally  dominate  the  dynamic  solution.  It  was  assumed  for  the  formulation  in 
SOM-TA  that  the  damping  ratio,  r„  was  constant  for  all  deformation  modes  of 
interest.  This  assumption  greatly  simplifies  the  solution  of  equation  (41) 
because  the  damping  coefficient  is  not  dependent  on  the  system  mass.  Equation 
(41)  can  then  be  simplified  to 


C  *  2K a  (42) 

where  the  stiffness  K  is  the  gradient  of  the  cushion  force-deflection  curve, 
and  a  is  a  constant  for  the  system. 

The  procedure  used  in  SOM-TA  to  calculate  the  current  value  of  the  damping  co¬ 
efficient  is  based  on  the  curent  value  of  K,  and  the  constant,  a,  which  is  de¬ 
termined  from  input  data.  The  slope  of  the  exponential  load/deflection  curve 

at  deflection  6  is 
c 

b<Sc 

K  =  Abe  c  (43) 

The  user  supplies  a  damping  coefficient,  C  ,  for  the  zero-deflection  condition 
which  the  program  uses  to  calculate  a  by  applying  equations  (43)  and  (42), 
thereby  resulting  in 

a  -  C° 

^  (44) 
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The  constant,  a,  and  the  current  gradient  of  the  cushion  load-deflection  curve, 
K,  are  used  at  each  time  step  to  determine  the  instantaneous  cushion  damping 
coefficient  from  equation  (42). 

2.4.2  Restraint  System  Forces 

The  method  used  in  calculating  the  forces  exerted  on  the  body  by  the  restraint 
system  differs  considerably  from  that  described  in  the  preceding  section  for 
the  contact  forces.  The  primary  reason  for  this  difference  is  that  the  re¬ 
straint  forces  do  not  act  at  any  fixed  points  on  the  occupant,  but,  rather, 
the  points  of  application  vary  wi  th  the  restraint  system  geometry. 

Although  other  configurations  can  be  selected  by  the  user,  a  restraint  system 
consisting  of  a  lap  belt  and  diagonal  shoulder  strap  will  be  used  as  an  example 
The  restraint  loads  are  transmitted  to  the  occupant  model  through  ellipsoidal 
surfaces  fixed  to  the  upper  and  lower  torso  segments.  These  surfaces  are  shown 
in  figure  13.  The  locations  of  the  anchor  points  A,,,  and  are  determined 
by  user  input  along  with  the  webbing  properties.  The  buckle  B  for  a  single 
shoulder  belt  is  located  according  to  an  input  parameter  which  specifies  the 
distance  from  the  appropriate  point,  in  this  case  Aj,  along  the  path  of  the 
lap  belt.  For  a  double-strap  shoulder  harness,  the  buckle  is  placed  on  the 
abdominal  contact  surface  between  its  intersections  wi th  the  thigh  surfaces. 
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Figure  13.  Restraint  system  configuration  variables. 


a2  =  R2 
b2  -  4/2 

c2  s  L2/2 

and  these  body  dimensions  are  defined  in  section  2.5, 

The  restraint  forces  are  determined  in  the  same  manner  for  both  the  upper  and 
lower  torso.  First,  the  belt  loads  are  calculated  from  the  displacements  of 
the  torso  segments,  and  the  resultant  force  on  each  segment,  is  then  applied  at 
the  point  along  the  arc  of  contact  between  the  belt  and  the  ellipsoidal  sur¬ 
face  where  the  force  is  normal  to  the  surface. 
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Explaining  this  procedure  in  further  detail  for  the  restraint  system  config¬ 
uration  shown  in  figure  13,  for  any  position  of  the  occupant,  the  coordinates 
of  the  left  shoulder,  the  hips,  and  the  buckle  connection  B  are  calculated  In 
the  aircraft  reference  frame.  The  length  of  each  side  of  the  lap  belt  is  equal 
to  the  sum  of  the  free  length  in  a  straight  line  from  an  anchor  point  (A^  or 
A2)  to  the  outermost  point  on  the  hip  contact  surface  (Cj  or  C2)  added  to  the 
length  of  an  arc  from  that  point  on  the  hip  to  the  buckle  (B).  The  shoulder 
belt  length  is  computed  as  the  sum  of  the  distance  from  the  anchor  point  (A3) 
to  a  point  of  tangency  on  the  top  of  the  shoulder  contact  surface  (D),  the 
distance  from  the  buckle  to  the  extreme  anterior  (forward)  point  on  the  ellip¬ 
soidal  chest  surface  (E),  and  the  length  of  an  arc  over  the  chest  surface  be¬ 
tween  points  D  and  E.  If  the  length  of  the  belt  segment  should  exceed  the 
equilibrium  (zero  load)  length  calculated  initially,  then  there  is  some  ten¬ 
sile  force  in  the  belt.  The  resultant  force  on  each  segment  is  the  vector  sum 
of  the  belt  forces.  Friction  between  the  shoulder  belt  and  chest  along  the 
length  of  the  belt  is  taken  into  account  by  reducing  the  load  in  the  belt 
between  the  chest  and  buckle  by  a  constant  fraction  of  the  load  in  the  free 
length  between  the  anchor  point  and  the  body  surface.  The  resultant  force  on 
the  lower  or  upper  torso  segment  may  be  written  generally  as 

I  =  F  xi  +  Fyj  +  F  zk  (47) 

where  F  ,  F  ,  and  F  are  components  in  the  local,  segment-fixed  coordinate 

A  J  Z 

system. 

To  find  the  point  on  the  segment  where  £_  is  normal  to  the  surface,  consider 
first  the  equation  of  an  ellipsoid: 

x2/a2  +  y2/b2  +  z2/c2  =  1  (48) 

which  may  also  be  expressed  in  functional  form  as 

f  (x,  y,  z)  =  x2/a2  +  y2/b2  +  z2/c2  -  1  (49) 

where  the  ellipsoid  can  be  regarded  as  the  level  surface  f=0  of  the  function. 
At  any  point  (x,  y,  z)  on  the  surface,  the  gradient  of  f  is  normal  to  the  sur¬ 
face.  The  gradient  is  given  by 
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grad  f  =  (2x/ac)  i  +  (2y/b‘)  J  +  (2z/efc)fe 
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and  at  the  point  of  application  of  the  resultant  force,  grad  f  is  colllnear 
with  £.  Making  use  of  the  proportionality  between  the  components  of  the  two 
vectors, 


i  T 


Fx  -  Cx/a 


Fy  -  Cy/b‘ 
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Fz  -  Cz/c‘ 


where  C  is  an  arbitrary  constant.  Solving  equation  (51)  for  the  coordinates 
(x,  y,  z)  and  suostituting  into  equation  (48) 


which  leads  to 


/F  a\z  /F  b\2  /F  c' 

/  x  .  y_\  .  z 
[nr  +  ^rj  +  pr 


r2  p  2  2  ,  r  2.  2  ,  r  2  2 

C  =  Fx  a  +  Fy  b  +  Fz  c 


if 

r* 


C  =  +  F  2b2  +  F  2c2 
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the  point  of  application  of  F_  is  then 


x  =  Fxa7C 


y  =  FybVC 


z  =  Fzc7C 


(55) 


wi  th 

C  «  -Jf  2a2  +  F  2b2  +  F  2c2 
f  X  y  z 

The  negative  sign  on  C  can  be  explained  by  the  fact  that  each  coordinate  in 
the  local  system  is  opposite  in  sign  to  the  corresponding  component  of  the  re¬ 
sultant  force,  or 

x  >  0  1  f  Fv  <0 

y  >  0  if  Fy  <  0  (56) 

z  >  0  if  Fz  <  0 

The  capability  of  the  belt's  point  of  application  of  resultant  belt  loads  to 
move  relative  to  the  torso  surfaces  allows  simulation  of  the  "submarining" 
under  the  lap  belt. 

2.5  OCCUPANT  DIMENSIONS  AND  INERTIAL  PROPERTIES 


Characteristics  required  by  the  occupant  model  for  each  of  the  segments  are 
the  length,  mass,  center  of  mass  location,  and  moments  of  inertia.  Also  re¬ 
quired  are  the  axial  and  flexural  stiffnesses  and  damping  coefficients  for 
both  vertebral  elements,  as  well  as  compliance  characteristics  for  the  chest 
and  abdomen.  It  is  assumed  that,  for  each  segment,  a  line  connecting  the 
joints  is  a  principal  axis,  so  that  the  required  moments  of  inertia  are  all 
principal  moments.  For  each  torso  element,  of  the  two-dimensional  model,  the 
center  of  mass  may,  in  general,  be  offset  from  the  spine  as  shown  in  figure  6. 
Moments  of  inertia  are  then  moments  with  respect  to  axes  located  at  the  mass 
center.  Final  data  required  to  describe  the  occupant  are  radii  of  26  contact 
surfaces,  which  are  ellipsoids  and  spheres. 

For  two  "standard"  occupants,  a  50th-percentile  human  male  and  a  50th-percentile 
anthropomorphic  (Part  572)  dummy,  all  the  required  data  are  stored  within  the 
program.  For  other  nonstandard  occupants,  the  above-described  data  must  be 
provided  as  input.  The  human  data  are  based  on  the  U.S.  Air  Force  drawing 
board  manikin  (reference  7),  whose  dimensions  were  extracted  by  the  FAA  Civil 
Aeromedical  Institute  and  tabulated  with  appropriate  inertial  properties  in 
reference  8.  Dummy  dimensions  and  characteristics  were  obtained  from  refer¬ 
ence  9. 
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Most  of  the  dimensions  and  inertial  properties  used  in  SOM-TA  were  taken  from 
references  8  and  9,  but  because  fewer  segments  are  used  in  the  mathematical 
model,  some  properties  needed  to  be  combined.  That  is,  In  the  model,  the 
hands  are  combined  with  the  lower  arms,  the  feet  with  the  lower  legs,  and  the 
mid  torso  with  the  lower  torso.  The  distance  from  the  wrist  pivot  to  the  mass 
center  of  the  hand  was  added  to  the  lower  arm  length;  a  contact  sphere  cen¬ 
tered  on  the  end  of  this  link  accounts  for  the  remaining  reach  of  the  hand. 
Similarly,  the  distance  from  the  ankle  pivot  to  the  foot  center  of  mass  was 
added  to  the  lower  leg  length,  and  a  contact  sphere  was  sized  to  provide  the 
correct  distance  between  the  knee  pivot  and  the  floor.  Composite  mass  center 
locations  and  moments  of  inertia  were  calculated  for  the  lower  torso,  lower 
arms,  and  lower  legs. 

2.5.1  Body  Segment  Dimensions 

The  basic  dimensions  of  the  occupant  segments  that  are  required  in  writing  the 
equations  of  motion  are  illustrated  in  figure  14.  The  lengths  of  the  segments 
are,  in  most  cases.,  effective  "link  lengths"  between  joint  centers,  rather  than 
standard  anthropometric  dimensions  based  on  external  measurements.  These 
lengths,  for  the  standard  50th-percent11e  occupants,  are  presented  in  table  1. 

For  segments  other  than  the  torso  segment,  the  distance  of  the  mass  center  of 
segment  n  from  the  end  nearest  the  body  reference  point  (m^)  i s  p n .  The  dis¬ 
tance  between  the  mass  center  and  the  far  end  is  given  by 
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Pn  =  Ln 


(57) 


Note  that  the  lengths  of  the  torso  segments,  and  L^,  are  not  used  in  the 
two-dimensional  model  shown  in  figure  5  although  they  are  used  to  generate 
contact  surface  ellipsoids  for  the  graphics  display.  However,  the  center  of 
mass  distances,  and  p2,  and  the  spinal  length,  S,  are  used  and 


L1  +  L2  =  P1  +  p2  +  S 


(58) 


Also,  the  seated  height  is  equal  to 


R1  +  pl  +  s  +  P2  +•  N  +  L3/2 
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TABLE  1.  BODY  SEGMENT  LENGTHS  (IN.) 


Segment 

50th-Percentile 

Air  crewmember 
(reference  8) 

Part  572 
Dummy 

(reference  9) 

Lower  Torso,  L^ 

9.44 

10.5 

Upper  Torso,  L2 

13.1 

11.5 

Neck,  N 

5.10 

4.88^ 

Head,  L3 

8.50^ 

8.35^ 

Upper  Ann,  L^ 

11.6 

11.3 

Lower  Arm,  L^ 

14.8^ 

1 3. 3^ 3  ^ 

Upper  Leg,  lQ 

17.1 

16.5 

Lower  Leg,  Lg 

18.4(3) 

18.0(3) 

Spine,  S 

12.4(3) 

10.85(3) 

Seated  Height 

3  7 . 0  ^ 3  ^ 

36.0(3) 

(1)  Scaled  from  manikin  drawing. 

(2)  Scaled  from  Part  572  drawing. 

(3)  Calculated. 


where  is  the  radius  of  the  pelvic  contact  ellipsoid  in  the  midsagittal  plane 
N  is  the  neck  length,  and  L^  is  the  head  segment  length  (z-di rection,  approxi¬ 
mately  menton  to  top  of  head,  in  anthropometric  terms). 

2.5.2  Body  Segment  Weights  and  Center  of  Mass  Locations 

!  ’  . -  —  ~  — . — . 

I 

Body  segment  weights  and  axial  locations  of  centers  of  mass  for  the  two  stan¬ 
dard  occupants  are  presented  in  table  2.  The  dimensions  e^  and  e^,  by  which 
the  torso  mass  centers  are  offset  from  the  spine  were  determined  to  be  0.20 
and  0.70  in.,  respectively,  based  on  simulation  of  vertical  (+G  )  tests.  A 
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I 
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H 

Livir. 


value  of  1.17  in.  for  the  corresponding  dimension  for  the  head,  e3,  Is  based 
on  Part  572  dummy  drawings. 

2.5.3  Body  Segment  Moments  of  Inertia 

Body  segment  moments  of  inertia  are  presented  in  table  3.  As  indicated,  mo¬ 
ments  of  Inertia  with  respect  to  lateral  (y)  axes  were  taken  from  references  8 
and  9  for  the  standard  human  and  dumir\y  occupants,  respectively.  The  moments 
of  inertia  with  respect  to  the  segment  x-  and  z-axes  were  determined  using  ap¬ 
proximations  to  segment  geometry.  The  torso  and  head  segments  were  approxima¬ 
ted  by  ellipsoids.  Assigning  appropriate  anthropometric  dimensions  to  the  el¬ 
lipsoid  axes,  ratios  Ix/Iy  and  I2/Iy  were  calculated  for  unit  mass.  These 
ratios,  multiplied  by  the  I  from  reference  8  or  9,  gave  values  of  Ix  and  Iz 
for  the  torso  and  head  segments.  The  Identical  procedure  was  used  for  the  ex¬ 
tremities,  except  that  these  segments  were  approximated  by  solid  circular  cyl¬ 
inders. 


TABLE  3.  BODY 

SEGMENT  MOMENTS  OF  INERTIA  (lb-in. -sec2) 

Seqment 

‘x 

. .  !y  . 

50th- Pe rcentile 
Alrcrewmember 
(reference  8) 

Part  572 

Dummy 

(reference  9) 

’z 

Lower  Torso 

4.03 

1.19 

0.760 

2.323 

Upper  Torso 

2.37 

3.29 

0.926 

1.70 

Neck 

- 

0.019 

0.0177 

- 

Head 

0.160 

0.199 

0.266 

0.233 

Upper  Arm 

0.131 

0.120 

0.135 

0.022 

Lower  Arm 
(including 

hand) 

0.105 

0.254 

0.185 

0.195 

Upper  Leg 

0.212 

1.41 

1.22 

0.873 

Lower  Leg 
(includi nq 

foot) 

1.28 

1.17 

0.994 

0.505 
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2.5.4  Body  Contact  Surfaces 

Twenty-six  surfaces  are  defined  on  the  body  for  calculation  of  external  forces 
exerted  on  the  occupant  by  the  seat  cushions,  the  restraint  system,  and  the 
rear  of  the  seat  in  front.  These  surfaces  are  ellipsoids,  cylinders  and  spheres, 
as  shown  in  figure  15,  The  dimensions  of  these  surfaces,  listed  in  table  4, 
were  obtained  from  anthropometric  data  in  references  8  and  9  or  scaled  off  the 
drawings  of  the  manikin  and  dummy.  The  surfaces  and  the  dimensions  required  for 
their  description  are  illustrated  in  detail  in  figure  16. 


Surface  Identification 

1.  Lower  Torso 

2.  Upper  Torso 

3.  Head 

4.  Right  Upper  Arm 

5.  Right  Forearm 

6.  Left  Upper  Arm 

7.  Left  Forearm 

8.  Right  Thigh 

9.  Right  Lower  Leg 

10.  Left  Thigh 

11.  Left  Lower  Leg 

12.  Right  Knee 

13.  Left  Knee 

14 .  Right  Foot 

15.  Left  Foot 

16.  Right  Hip 

17.  Left  Hip 

18.  Right  Shoulder  Joint 

19.  Left  Shoulder  Joint 

20.  Right  Elbow 

21.  Left  Elbow 

22 .  Right  Hand 

23.  Left  Hand 

24.  Neck 

25.  Right  Shoulder 

26.  Left  Shoulder 


Figure  15.  Occupant  contact  surfaces 


TABLE  4.  BODY  CONTACT  SURFACE  RADII  (IN.) 

Contact 

Surface 

50th- Percentile 

Ai rcrewmember 
(scaled  from  manikin) 

Part  572 

Dummy 

(scaled  from  drawinq) 

Lower  Torso 

4.00 

4.50 

Upper  Torso 

5.00 

4.50 

Neck(1) 

2.00 

2.00 

Head^ 

3.75 

3.44 

Upper  Arm 

2.10 

1.95 

Lower  Arm 

1.65 

1.85 

Upper  Leg 

3.20 

3.10 

Lower  Leg 

2.23 

2.30 

Hip^ 

3.56 

3.56 

Shoulder^ 

2.00 

2.00 

Foot 

1.60 

1.60 

(1)  Neck  circumference  divided  by  2 tt  . 

(2)  Head  length  (anterior -posterior ) . 

(3)  Hip  breadth  ( s 1 tti ng )/ 2  -  hip  link  length  (Ln). 

(4)  Shoulder  breadth/2  -  shoulder  link  length  (L^). 


2.5.5  Joint  Rotation 

The  results  of  several  studies  on  the  limits  of  human  joint,  motion  have  been 
published.  Two  of  these  studies,  in  particular,  were  examined  for  applicabil¬ 
ity  to  the  occupant  model.  First  of  all,  Dempster's  (reference  10)  data  on 
link  lengths  and  inertial  properties  were  used,  as  discussed  in  preceding  sec¬ 
tions,  so  it  was  considered  appropriate  to  include  his  joint  data  here.  Gian- 
ville  and  Kreezer  (reference  11)  presented  limits  of  joint  motion  for  both 
voluntary  and  forces  rotation;  their  results  appear,  along  with  Dempster's,  in 
table  5.  Definitions  of  the  various  joint  motions  are  illustrated  in  figure  17. 
Also  included  in  table  5  are  the  rotations  required  for  the  Part  572  anthropo¬ 
morphic  dummy. 
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♦Transverse  plane 


Figure  17.  Motion  diagrams. 


All  of  the  rotations  possible  in  the  mathematical  model  are  included  in  table  5 
and  figure  17,  but  some  are,  naturally,  more  important  than  others  in  determin¬ 
ing  permissible  ranges  of  motion  for  the  me ‘el .  For  the  head,  ventriflexi on 
(B)  is  certainly  the  most  important  component  of  motion  for  frontal  impact. 
Dorsi flexion  (A)  may  also  be  important  for  frontal  impact,  but  the  angles  re¬ 
ported  are  sufficiently  chose  to  those  for  ventri  flexion  to  be  considered  the 


same.  Laterial  flexion  (C)  is  certainly  less  important  since  a  pure  lateral 
Impact  of  an  aircraft  would  be  rare  indeed,  and  rotation  (D)  will  have  an  in¬ 
significant  effect  on  model  response.  Therefore,  the  limiting  rotation  Be 

bi 

(see  section  2.3)  for  the  neck  joint  (i  *  2)  has  been  taken  as  the  limit  for 

voluntary  ventriflexion,  or  60  degrees.  This  angle  can  be  found  In  table  6, 

along  with  the  limiting  angles  for  the  other  body  joints.  For  all  of  the  other 

angles,  flexion  is  the  most  significant  component  for  the  type  of  motion  that 

can  usually  be  expected  to  take  place  in  a  crash  environment.  Therefore,  the 

limiting  angles  were  all  taken  as  the  limits  for  voluntary  flexion.  Note  that, 

for  the  hip  joint,  the  reference  position  of  the  body  used  in  the  mathematical 

model  Includes  90-degree  flexion.  Therefore,  this  amount  has  been  subtracted 

from  the  angle  reported  in  table  5,  which  is  defined  relative  to  the  standard 

anatomical  reference  position.  Since  the  seated  position  appears  to  aid  in 

flexion  of  the  hip  joint,  the  largest  angle  in  the  table,  the  one  given  by  the 

Part  572  specifications,  was  used  in  determining  3$  which  is  thus  given  by 

Br  =  120°  -  90°  =  30°.  7 

/ 


TABLE 

6.  JOINT  LIMITING  ANGLES 

Joint 

Location 

Angle  -  Bo 
(deg)  i 

1 

Back 

40 

2 

Neck 

60 

3,  5 

Shoul der 

180 

4,  6 

El  bow 

142 

7,  9 

Hi  p 

30 

8,  10 

Knee 

125 

2.5,6  Body  Stiffnesses 

The  lumbar  spine  and  neck  of  the  two-dimensional  model  possess  exponential 
stiffness  characteristics,  in  the  form  of  equation  (40),  for  both  axial  and 
rotational  deformation.  Exponential  stiffness  characteristics  for  the  abdo¬ 
men  and  chest  are  used  to  soften  the  input  force-deflection  values  for  the 
lap  belt  and  shoulder  belt,  respectively. 


3.0  SEAT  MODEL 


The  seat  structure  Is  modeled  using  the  finite  element  method  of  analysis. 

This  method  has  been  selected  because  It  Is  not  reliant  on  previous  testing, 
and  it  has  the  flexibility  to  deal  with  a  wide  range  of  design  concepts.  The 
specific  finite  element  formulation  used  in  the  program  is  based  on  the  WRECKER 
II  program,  developed  at  the  ITT  Research  Institute  (reference  12). 

Program  SOM-TA  seat  analysis  uses  three-dimensional  beam  elements.  It  has  the 
capability  to  model  large  displacements,  nonlinear  material  behavior,  local 
buckling,  and  various  end  release  conditions  for  beam  elements.  The  large 
displacement  formulation  separates  the  element  displacement  field  Into  a  rigid 
body  rotation  and  translation  associated  with  a  local  coordinate  system  that 
moves  with  the  element  and  small  element  distortions  relative  to  the  current 
position  of  the  element  coordinate  system.  This  formulation  can  accommodate 
extremely  large  rotations  and  deflections  with  accuracy  depending  on  the  size 
of  the  elements  relative  to  the  curvature  of  the  structure.  Nonlinear  material 
formulation  is  based  on  a  uniaxial  elastic-plastic  stress-strain  law.  beam 
end  release  conditions  include  axial  load,  shear,  moment,  and  torque  releases. 
Also,  a  simple  local  buckling  model  for  thin-walled  tubes  subjected  to  axial 
compressive  and/or  bending  loads  was  incorporated  into  the  program.  This  model 
simulates  the  reduction  in  bending  rigidity  of  the  tube  as  the  cross  section 
distorts  during  local  buckling. 

3.1  SOLUTION  PROCEDURE 


The  solution  procedure  is  based  on  formulation  of  equations  of  quasi-static 
equilibrium  for  the  finite  element  model  in  an  incremental  form  according  to 


-T  (7j'i+l  -  V  =  Fi+1  "  F1 


(60) 


where 


=  Tangent  stiffness  matrix 

-  The  vector  cf  external  applied  loads  at  i+lth  solution  time 
step. 

T?  =  The  vector  of  internal  forces  associated  with  element  deforma- 
1  th 

tions  at  1  solution  time  step. 
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—  J.  L. 

u.  -  D1  s placements/ rotations  at  1  n  solution  time  step. 

The  external  forces,  F*\  including  restraint  system  loads,  occupant  loads  on 
the  seat  pan  and  seat  back,  and  seat  support  reactions  are  treated  as  static 
loads  on  the  seat  structure.  The  mass  of  the  seat  structure  is  neglected, 
since  in  most  simulations  it  will  be  a  small  percentage  of  the  total  occupant 

weight,  and  the  seat  structure  is  assumed  to  be  in  a  quasi-static  equilibrium 

.  — E 

with  the  applied  external  forces  F  . 

The  tangential  stiffness  matrix  depends  on  the  state  of  stress  of  the  seat 
structure  and  varies  with  time  during  the  simulation.  Therefore,  it  must  be 
recomputed,  assembled  and  inverted  at  selected  time  steps  during  the  solution. 
The  major  computational  effort  is  the  Inversion  of  which  requires  NB  multi¬ 
plications  where  3  is  the  serni bandwidth  and  N  is  the  number  of  degrees  of  free¬ 
dom. 


f 

i 

♦ 

i 
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3.2  BEAM  ELEMENT  FORMULATION 


Large  displacement  formulation  separates  the  beam  element  di spl acement  field 
into  rigid  body  rotation  and  translation  and  small  element  distortion  relative 
to  the  current  position  of  the  element  coordinate  system.  After  the  rigid  body 
motion  is  removed,  it  is  possible  to  use  the  classical  small  deformation  finite 
element  formulations.  Consequently,  extremely  large  rotations  and  translations 
can  be  accommodated  with  accuracy  depending  on  the  size  of  the  elements  relative 
to  the  deformed  curvature  of  the  structure. 

The  beam  element  is  based  on  the  conventional  small -deflection  formulation  in¬ 
volving  cubic  displacement  fields  for  transverse  displacements  and  linear  dis¬ 
placements  for  axial  and  torsional  displacements. 

From  the  principle  of  virtual  work,  the  general  form  of  the  beam  element  tan¬ 
gent  stiffness  matrix  (reference  12)  is  given  by 

Kt  =  /  DT  C  D  dV  (61) 

-T  v 

where  K.^  =  Element  tangent  stiffness  matrix 

l)  =  Matrix  that  relates  element  strains  to  nodal  dis¬ 
placements 

C  -  Constitutive  matrix  that  relates  stresses  and  strains. 
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For  linear  elastic  beams,  the  form  of  the  equation  given  in  equation  (61)  re¬ 
sults  In  the  classical  12  x  12  elastic  beam  stiffness  matrix  shown  in  fig¬ 
ure  18. 

Where  A  -  area  of  cross  section 

E  =  modulus  of  elasticity 
I  a  moment  of  inertia 
G  =  modulus  of  rigid!  ty 
J  =  torsion  constant 
L  =  length  of  member 

(Ny.)»  (M  )  °  moment  about  x,  (y),  (z)  at  node  1 

K  .  (Fv  )»  (F,  )  =  force  in  xt  (y),  (z)  direction  at 
xi  '1  Z1  node  1 

0x  *  (0v  <0z  )  =  rotat1’on  about  x,  (y),  (z)  axis  at 

l  1  node  i 

ux  .  (v  ),  (w  )  ■  x,  (y),  (z)  component  of  displacement 
i  yi  zi  at  node  i 


For  nonl  inear  materials  the  constitutive  matrix  C.  is  given  by 

r  Acr 
-  ~  Ac 

where  Aa  -  incremental  stress 
Ac  =  Incremental  strain. 


(62) 


Using  the  cubic  shape  functions  for  the  transverse  displacements  and  linear 
shape  functions  for  axial  and  torsional  displacements  in  matrix  D  It  can  be 
shown  that  (reference  12) 
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where  £  3  element  length 


*  6x  -  2£ 


=  6x  -  4£ 


;  |2.  dA 
A  Ae 


6V  =  /  y  4“  dA,  6  =  /  Z  ~  dA 

y  A  Ae  z  A  Ae 


.  -  Ao  ,  2  Aa 

«  =/yz  jj-dfl,  6  =  {y  -££  dA, 


j  r  2  Aa 

sz  -  '  2  li  dA 


The  Integrals  defined  in  equations  (63)  and  (64)  are  calculated  numerically 
through  the  cross  section  and  along  the  length  of  the  beam  element.  The  inte 
gration  is  piecewise  linear  through  the  depth  and  linear  along  the  length. 


,w*r 


Once  the  element  stiffness  matrices  have  been  evaluated,  it  is  necessary  to 
assemble  each  matrix  to  some  global  system.  To  accomplish  this,  an  element 

A  A  ^ 

coordinate  system  (x,  y,  z)  is  attached  to  each  beam  element  and  serves  to 
define  the  rigid  body  rotation  and  translation  of  the  element.  The  orien¬ 
tation  of  element  axes  (x,  y,  z)  with  respect  to  the  global  axes  at  any  time 
is  established  by  the  components  of  three  vectors  e.,  e_,  e-  which  remain  fixed 

AAA  1  U 

along  the  element  axes  (x,  y,  z),  respectively,  as  the  element  translates  and 
rotates.  If  these  three  unit  vectors  e^,  e2>  e^  form  the  columns  of  a  3  x  3 
matrix  _£,  then  any  vector  v  can  be  transformed  from  element  to  global  coordinate 
system  by  the  following  time  dependent  transformation: 

vg-£ve  <65> 

For  beam  elements  (figure  19),  the  x-axis  is  defined  by  a  line  connecting  the 
end  points  of  the  beam,  and  the  y-axis,  by  a  line  normal  to  the  x-axis  and 
lying  in  a  plane  containing  both  the  coordinate  reference  point  and  x-axis. 

The  remaining  z-axi s  is  determined  as  a  normal  to  x-  and  y-axes  by  the 
righthand  rule. 

Also  accounted  for  in  the  transformation  between  the  local  and  global  systems 
is  the  offset  of  the  beam  end  nodal  point  from  the  beam  cross-section  centroid 
and  shear  center  (figure  19b).  This  transformation  Is  particularly  important 
for  open  beam  cross-sections  where  the  centroid  and  the  shear  center  do  not 
coincide.  The  appropriate  transformation  for  the  beam  end,  consistent  with 
the  beam  kinematic  assumptions,  can  be  written  as: 


10  0 


vi  >  -  0  10 


wj  0  0  1 


0  -(z  -  z  )  +(y  -  y  ) 

n  eg'  eg 


-»n  0 


l 


a 


where  u.. ,  v^,  =  Displacements  at  the  shear  center  In  beam  local  x,  y, 

and  z  directions. 


u  ,  v  ,  w  «  Displacements  at  the  node  point  in  beam  local  x,  y, 
n  ri  n 


i  1  1  and  z  directions. 


A 


a 


6,6,0  =  Cross-section  rotations  about  beam  local  x,  y,  and  z 

x  y  z 

1  i  i  axes. 


y  ,  z  =  Distance  between  the  shear  center  and  the  node  point 
n  n 


(see  figure  19b) . 


fA 


y  ,  Z  =  Distance  between  the  shear  center  and  the  centroid. 
Jcg’  eg 


3.3  INTERNAL  RELEASES  IN  BEAM  ELEMENTS 


The  program  allows  the  use  of  "end  releases"  in  beam  elements.  The  beam  end 
release  conditions  include  axial  load,  shear,  moment,  and  torque  releases 
at  the  end  nodes.  Combinations  of  these  releases  are  possible,  e.g.,  a  ball 
joint  is  formed  by  releasing  moments  in  three  different  directions. 

Depending  on  the  direction  of  the  end  release,  the  original  beam  element  stiff 
ness  matrix  is  modified  as  follows: 

h  \  ■  f'e  <6?> 

where  K  -  Element  stiffness  matrix 
— e 

^  =  Element  end  displacements/rotations 

F  ■  Element  end  forces/moments 
—e 

Rewriting  equation  (67)  as 


where  Ap  =  Fixed  end  displacements 

Ar  =  Released  end  displacements 
Fp  =  Fixed  end  forces/moments 
From  equation  (68) 

J<FF  +  ^FR  ^F 


(68) 


The  modified  beam  element  stiffness  matrix  K  *  is 

— € 


then  given  by: 


-FF  '  ^FR  — RR  K-RF 


0 

0 


(71) 


3.4  MATERIAL  NONLINEARITIES 

Nonlinear  material  formulation  for  beam  elements  is  based  on  uniaxial  elastic 
plastic  stress-strain  law. 


The  computational  algorithm  is  based  on  the  work  of  Hartzman  and  Hutchinson, 
(reference  14),  as  specialized  for  small  strain  and  plane  stress  conditions. 

The  current  stress  state  at  a  point  on  the  cross  section  of  the  beam  is  estab¬ 
lished  as  follows:  let  the  prior  stress  state  of  the  point  under  consideration 
be  o^O  and  a  small  increment  in  strain  from  the  prior  state  to  the  current 
state  by  Ac  .  First  a  tentative,  current  stress  state,  o  is  calculated  as 

X  X 

though  the  strain  increment  were  completely  elastic. 


(72) 


where  E  =  Modulus  of  Elasticity 
Y  =  Pol sson' s  ratio 

An  effective  stress  o  =  o  is  used  to  determine  whether  plastic  flow  has  taken 

ex 

place  during  the  strain  increment.  If  a  is  less  than  the  prior  effective 
G  ^ 

stress  ag  at  which  yielding  occurred,  elemental  loads  are  decreasing  and  the 

tentative  stress  calculated  a  is  the  correct  value  of  the  current  stress  o  . 

x  0  x 
If  o  is  greater  than  the  prior  effective  stress  o  ,  at  which  yielding  oc- 

curred,  tentative  stresses  calculated  must  be  modified  to  account  for  the 

plastic  behavior.  Hartzman  and  Hutchinson  (reference  14)  have  shown  that  the 

true  value  of  the  current  effective  stress  will  be: 


o/ 


The  current  state  of  stress  Is  then  given  by 


a 


-  JlLA  n 
x  “  r+'S'A  x 


(74) 


where 


-1) 


3.5  LOCAL  BUCKLING  OF  BEAM  ELEMENTS 

Local  buckling  Is  one  of  the  failure  modes  of  thin-walled  tubes,  typically  used 
in  light  aircraft  seats,  when  subjected  to  axial  compressive  and/or  bending 
loads.  Evaluation  of  the  results  of  dynamic  tests  of  seat-occupant  systems 
have  indicated  that  the  local  buckling  of  thin-walled  tubes  have  a  significant 
effect  on  the  response  of  the  seat-occupant  system.  Therefore,  a  simple  local 
buckling  model  for  thin- walled  tubes  was  incorporated  into  the  program. 

Local  buckling  Is  traditionally  expressed  in  terms  of  a  moment-bending  curva¬ 
ture  diagram  as  shown  in  figure  21. 


Curvature  change  -  9 


Buckled 


Figure  21.  Moment  capability  versus  curvature 
of  a  thin-walled  circular  tube. 


The  cross  section  goes  through  several  stages  of  deformation  as  the  structure 
bends.  Although  during  this  cross-section  distortion  axial  stresses  redistrib¬ 
ute  themselves.  It  was  proposed  In  reference  15  that  the  reduced  bending  rigid¬ 
ity  is  most  strongly  related  to  the  loss  of  lateral  moment  arm  of  the  axial 
fo rces. 

For  a  circular  tube  under  uniform  compression  the  stress,  at  which  local 
buckling  starts  is  given  in  reference  16  as 

=>L'KcE^>  ™ 


where 


0.4  <  K  <  1.2  (recommended  values) 


E  =  modulus  of  elasticity 
t  =  wall  thickness  of  the  tube 
D  =  diameter  of  the  tube. 


It  was  also  proposed  in  reference  16  that  for  tubes  subjected  to  bending,  the 
local  buckling  stress,  oL>  as  a  function  of  material  and  geometry,  can  be  taken 
as  that  for  a  circular  tube  under  uniform  axial  compression,  since  in  bending, 
a  significant  portion  of  the  circumference  is  subjected  to  a  relatively  uniform 
compression  field.  However,  dynamic  tests  of  seat-occupant  systems  have  indi¬ 
cated  that  the  thin-walled  circular  tubes  have  suffered  local  buckling  at  stress 
levels  much  below  that  predicted  by  equation  (75).  This  was  partly  due  to  the 
fact  that  the  bending  stresses  have  exceeded  the  yield  point  and  the  tubes  have 
also  suffered  plastic  deformation. 


To  account  for  plasticity  as  well  as  cross-sectional  properties,  an  empirical 
relation,  equation  (76)  was  used  in  he  program  to  predict  the  local-buckling 
stress 


0.4  E  (t/D) 


0  ^  -  Minimum 


(76) 
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where  =  local  buckling  stress 

o  3  compressive  yield  stress 
t  3  wall  thickness  of  the  tube 
D  =  diameter  of  the  tube.* 

The  cross  sections  of  beam  elements  are  defined  as  thin-walled  plate  segments 
by  specifying  their  end  points  and  thickness  as  part  of  the  input  data  (fig¬ 
ure  22).  The  tangent  stiffness  matrix  for  nonlinear  beam  elements  is  then  com¬ 
puted  by  evaluating  the  Integrals  given  in  equation  (66)  over  the  plate  segments 
defining  the  beam  cross  section. 

The  program  computes  the  stresses  at  all  segment  end  points  across  the  cross 
section  at  each  end  of  the  beam  elements.  These  stresses  are  then  compared 
with  the  local  buckling  stress,  o^,  computed  from  equation  (76)  for  each  beam 
element.  If  any  of  the  compressive  stresses  at  the  segment  end  points  exceed 
o^,  the  deformation  of  the  cross  section  is  modeled  by  modifying  the  radial 
location  of  that  segment  end  point  using  the  following  expression 


z 


Figure  22.  Circular  tube  cross  section  defined 
by  eight  thin- walled  plate  segments. 


*For  rectangular  tubes  the  diameter  of  a  ci rcumscribi ng  circular  tube  is  used. 


1 


'll 


tfhere 


radial  location  of  the  Gaussian  station  at  time 
step  1 


aL  ■  local  buckling  stress 

ai+1  »  compressive  stress  at  the  segment  end  point  at  time  step  1+1 
K  ■  local  buckling  constant  (0.50  recommended). 


Consequently,  reduced  bending  rigidity  of  the  cross  section  due  to  the  decrease 
in  the  lateral  moment  arm  of  axial  forces  during  local  buckling  can  be  modeled. 


Iff’.' 
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4.0  SIMULATION  COMPUTER  PROGRAM 


fa 


The  digital  computer  program  based  on  the  occupant  and  seat  models  described 
in  sections  2  and  3  is  called  Seat/Occupant  Model  -  Transport  Aircraft  (SOM-TA) 
It  has  been  written  entirely  in  FORTRAN  to  ensure  a  high  degree  of  compatibil¬ 
ity  with  various  digital  computer  systems.  During  development,  the  program 
has  been  run  on  IBM  and  Digital  computer  systems. 

The  elements  of  the  program  can  be  considered  in  terms  of  three  general  opera¬ 
tions: 

•  Input  and  initialization 

•  Solution 

•  Output 

which  are  summarized  below  and  discussed  in  detail  in  the  sections  following. 
The  general  flow  of  the  program  is  illustrated  in  figure  23.  Input  data  de¬ 
scribing  the  occupant  and  crash  conditions  are  read  first.  Finally,  the  seat 
data,  either  simple  dimensions  describing  a  rigid  seat  model  or  detailed  de¬ 
sign  data  for  the  finite  element  seat  analysis,  are  provided.  Based  on  the 
input  data,  the  values  of  constants,  such  as  occupant  dimensions  and  proper¬ 
ties  are  calculated,  and  the  initial  position  of  the  occupant  is  determined. 


K 


The  solution  loop  is  entered  for  the  first  time  with  the  aircraft  initial 
velocity  and  the  occupant  initial  position.  At  each  subsequent  entrance  to 
the  loop,  the  current  aircraft  displacement,  velocity,  and  acceleration  com¬ 
ponents  are  calculated.  The  equations  of  motion  for  the  occupant  are  set  up 
and  solved.  If  a  finite  element  seat  model  is  being  used,  the  forces  applied 
to  the  seat,  such  as  the  cushion  forces,  are  provided  to  the  seat  routines  for 
computation  of  seat  displacements.  At  time  increments  equal  to  a  predeter¬ 
mined  print  interval,  the  output  variables  requested  by  the  user  are  stored 
for  printing  after  completion  of  the  solution. 
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4.1  PROGRAM  INPUT 


Input  data  are  read  by  the  program  in  the  following  six  blocks: 

1.  Simulation  and  output  control  information 

2.  Cushion  properties 

3.  Restraint  system  description 

4.  Crash  conditions 

5.  Occupant  description 

6.  Seat  design  information. 


4.1.1  Simulation  Control  Information 

The  first  block  of  data  contains  the  information  required  for  controlling 
execution  of  the  program.  The  selection  of  system  of  units  (SI  or  English), 
two-  or  three-dimensional  occupant,  the  type  of  seat  (single-,  double-,  or 
triple-occupant)  along  with  identification  of  the  positions  that  are  occupied, 
secondary  Impact  Information  (if  desired),  total  simulation  time,  and  identi¬ 
fication  of  desired  output,  are  provided  here. 


4.1.2  Cushion  Description 


Cushion  load-deflection  characteristics  are  described  by  an  exponential  func¬ 
tion,  whose  coefficients  are  provided  as  input  data.  The  equilibrium  (zero 
load)  thickness  for  both  the  seat  and  back  cushions  are  also  given.  The  cush¬ 
ion  damping  coefficient  for  zero  deflection  described  in  section  2.4.1  is  also 
entered. 


4.1.3  Restraint  System  Description 


The  restraint  system  used  in  the  simulation  may  consist  of  a  lap  belt  alone  or 
combined  with  a  single-  or  double-strap  shoulder  harness.  The  webbing  force- 
elongation  curve  is  approximated  by  three  linear  segments,  which  are  described 


by  input  of  points  on  the  curve.  The  force  is  computed  by  linear  interpola¬ 
tion  In  this  table,  as  described  in  section  4.3.1.  The  slack  in  the  webbing 
is  also  provided  by  input  in  units  of  length. 

The  anchor  points  for  the  lap  belt  and  shoulder  harness  are  located  by  input 
of  rectangular  coordinates  in  the  aircraft  reference  system.  For  a  double¬ 
strap  shoulder  harness,  the  buckle,  or  point  of  connection  to  the  lap  belt,  is 
assumed  located  on  the  mid-point  of  the  lap  belt.  For  a  single  shoulder  belt, 
which  may  pass  over  either  the  left  or  the  riyht  shoulder,  an  input  parameter 
locates  the  buckle  by  the  length  of  webbing  between  the  buckle  and  the  lap  belt 
anchor  point.  This  length  may  be  zero  If  the  buckle  attaches  directly  to  a 
rigid  anchor  point. 

4.1.4  Crash  Conditions 

The  aircraft  crash  conditions  are  defined  by  the  initial  velocity  and  attitude 

and  the  acceleration  as  a  function  of  time.  Six  components  of  velocity  are 

required:  three  translational  in  the  aircraft  coordinate  system  (VY  ,  VY  ,  V7  ) 

*A  yA  Lk 

and  the  yaw,  pitch,  and  roll  rates  ($A,  eA.  <t>A).  Each  of  the  six  acceleration 
components,  which  define  the  acceleration  of  the  aircraft  coordinate  system, 
is  described  by  sixteen  points  in  time  and  acceleration.  An  example  of  an  ap¬ 
proximation  to  an  actual  acceleration  pulse  is  illustrated  in  figure  24.  Al¬ 
though  many  of  the  higher  frequency  oscillations  observed  in  the  actual  pulse 
probably  contribute  little  to  the  overall  response  of  the  occupant,  the  use  of 
a  large  number  of  points  reduces  the  effect  of  the  investigator's  subjectivity 
in  the  approximation, 

4.1.5  Occupant  Description 

Because  it  has  been  assumed  that  the  principal  user  of  this  program  is  inter¬ 
ested  primarily  in  the  seat  or  restraint  system,  a  minimum  of  information  is 
required  to  describe  the  occupants.  Data  for  standard  human  and  dummy  occu¬ 
pants,  as  described  in  section  2.5,  are  stored  within  the  program.  Additional 
data  must  be  provided  for  nonstandard  occupants. 
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TIME  -  SEC 


rigure  24.  Piecewise  approximation  to  aircraft 
acceleration  component. 


4.1.6  Seat  Design  Information 


The  input  data  required  to  describe  the  seat  consist  of  nodal  coordinates, 
material  properties,  cross-section  geometries,  element  locations,  arid  attach¬ 
ment  conditions,  as  described  earlier.  For  possible  use  in  restraint  system 
analyses  where  detailed  seat  response  may  not  be  important  or  seat  design 


4>1. 7  Seat  Back  Contact  Surfaces 

the  surfaces  Illustrate*  In  fi.»«  «■  ‘  a'  these  surfaces,  which 

to  specify  the  force-hefomation 

represent  cushion,  tray  table,  an*  amrest  surfaces. 


Figure  *'a 


-eat  back  contact  surfaces. 
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4.2  OCCUPANT  INITIAL  POSITION 


The  initial  position  of  the  aircraft  occupant  is  computed  from  the  input  pa¬ 
rameters  shown  in  figure  2b.  The  occupant  is  seated  at  Y-coordinate  of  the 
mid-plane  (plane  of  symmetry)  for  each  occupant  as  specified  in  the  input 
data.  The  angular  coordinates  Y..  (1  =  1,  2,  3,  4)  define  the  rotation  of 
segments  1-4  relative  to  the  ZA  axis  and,  because  of  the  symmetry  condition, 
segment  6  is  parallel  to  segment  4.  Positive  angles  are  shown  in  figure  26 
although  Y^  and  Y^  are  usually  negative,  considering  the  torso  to  be  approxi¬ 
mately  parallel  to  an  aftward-sl oping  seat  back.  The  angle  Yg  describes  the 
position  of  the  forearms  relative  to  the  upper  arms,  and  is  the  initial  value 
of  ag  and  a.,.  The  distance  X^  Is  the  initial  X-coordinate  of  the  heels  (the 
inferior  ends  of  segments  9  and  11).  The  procedure  described  below  consists 
of  seating  the  occupant  in  such  a  position  that  static  equilibrium  is  achieved 
among  the  forces  exerted  by  the  seat  cushion,  floor,  and  either  the  restraint 
systan  or  the  back  cushion. 

The  first  step  in  determining  the  initial  position  for  the  three-dimensional 
occupant  involves  calculating  the  Euler  angles  for  the  torso,  head,  and  arm 
segments,  since  this  procedure  does  not  require  consideration  of  the  forces 
due  to  the  cushions  and  floor.  Because  the  input  parameters  illustrated  in 
figure  26  define  the  position  of  the  occupant  in  the  aircraft  coordinate  sys¬ 
tem,  the  orientation  of  the  aircraft  must  be  described  in  the  inertial  system. 
For  an  aircraft  in  level  flight  with  cero  pitch,  roll,  and  yaw,  It  is  assumed 
that  the  aircraft  coordinate  axes  (X^,  Y^,  Z^)  are  parallel  to  the  fixed  co¬ 
ordinate  axes  (X,  Y,  Z.)  at  the  Initial  time  (t  =  0).  A  general  orientation  of 
the  aircraft  reference  frame  is  obtained  by  the  same  sequence  of  rotations 
defined  in  section  2.1.1  for  the  body  segments.  Defining  the  rotations 


Fi ciure  2.6,  Initial  position  input  parameters. 

t ne  orientation  of  the  aircraft  relative  to  the  inertial  system  is  described 
by  the  coordinate  transformation 
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(78) 


where  the  elements  of  [A]  are 


Ajj  *  cos  \1>A  cos  6A 

A12  «  cos  i|>A  sin  0A  sin  4>A  -  sin  il»A  cos  $A 

A13  =  cos  sin  0A  cos  $A  +  sin^A  sin  4>A 

A21  =  sin  ^A  cos  0 A 

A22  =  ^A  sin  0A  s’n  ^A  +  C0S^A  cos  <t>  a 

A23  =  sin  \PA  sin  0A  cos  4>A  -  cos  ^A  sin  <t>A 

A31  *  -s1n  6A 
A32  =  cos  0A  sin  4>a 

A33  =  cos  0  A  cos  <}>  A  (79) 

The  rotation  of  body  segment  n  relative  to  the  aircraft,  remembering  that  the 
symmetry  condition  requires  that  yn  is  parallel  to  YA,  is  described  by 


Combining  equations  (78)  and  (80)  results  in  the  angular  relationship  between 
the  local  coordinate  system  of  segment  n  and  the  inertial  system  expressed  by 
the  following  transformation,  which  is  a  function  of  the  input  Yn  and  the  air 
craft  pitch,  roll,  and  yaw: 


where  [B  j  is  given  by 


[Bn]  *  A 


cos  y  0  sin  y 
1  n  '  n 


0  1 


-sin  Tn  0  cos  “Y 


so  that  its  elements  are 


Bll*AllCOSYn  -A13sinYn 


B12  =  A12 


B13  b  A11  s^n  Y  n  +  A13  cos  y  n 
B"l  ’  A21  A53 


B  =  A 
22  rt22 


n  23  7n 


B23  =  A21  sln  Yn  +  A23  005  Yn 


831  ‘  A31  COSYn  •  A33  s1nYn 


B  =  A 
u32  m32 


B33  =  A31SinYntA33  COSYn 


Comparison  of  equation  (81)  with  equation  (1)  points  out  that  the  transforma¬ 
tion  matrices  [Tn]  and  [Bn]  are  equivalent.  Because  [Tn]  is  a  function  of 
Euler  angles  for  segment  n,  equating  the  elements  of  [Tn]  and  [Bn]  through 


ft; 


[Tn]  »  [Bn]  n  =  1,  2,  3,  4,  6,  8,  10  (84) 

permits  calculation  of  the  Initial  values  of  the  generalized  coordinates  from 
input  parameters  y  and  6^,  and  4> The  procedure  as  used  in  Program 
SOM-TA  is  outlined  below. 


First  0n  is  determined  as  follows: 


which  gives 


t"31  -  b"31  or  -Sine  „  =  B"31 


6  n  '  s*"'1 


The  roslne  is  then  found  by 

1  n 

cos  e  n  =  cos  [sin"  (-B^)] 
so  that  \pn  can  be  determined: 

T11  =  B11  or  cos  ^  n  cos  9  n  =  B11 


(85) 


gives 


»l»n  =  cos"1  (B^/cos  0n) 

and,  for  determination  of  <j>  , 

T1!-  =  B1!-  or  cos  0  cos  <j>  =  B-., 

33  33  n  n  33 


gives 


4> 


n 


cos 


-1 


(B33/cos  e n) 


(86) 


(87) 


Equations  (85)  through  (87)  ?.ce  used  for  segments  1,  2,  3,  and  4;  the  symmetry 
requirement  provides  the  Euler  angles  for  segment  6,  which  are  equal  to  those 


73 


for  segment  4.  At  this  point  the  generalized  coordinates  through  q2Q  have 
been  determined.  The  next  step  involves  seating  the  occupant  and  calculating 
Xj_ •  Y1 ,  and  Zj  (qj,  q2,  and  q3)  from  static  equl librium.'1' 

Because  the  problem  of  seating  the  occupant  is  statically  indeterminate,  cer¬ 
tain  simplifying  assumptions  are  made.  The  first  assumption,  which  is  approxi¬ 
mately  correct  for  typical  seating  positions,  is  that  15  percent  of  the  occu¬ 
pant's  weight  is  supported  by  the  floor.  In  other  words,  85  percent  is  sup¬ 
ported  by  the  seat  cushion  and,  depending  on  the  aircraft  attitude,  the  re¬ 
straint  system  or  the  back  cushion. 

A  first  approximation  to  the  initial  position  is  made  for  the  assumption  of 
level  flight  (@A  3  <1^  “  ^  3  O)*  The  cushion  forces  act  on  the  body  as  shown 
in  figure  27,  where  it  is  assumed  that  15  percent  of  the  occupant  weight  is 
supported  by  the  floor,  as  discussed  in  the  preceding  paragraph.  Summing 
forces  gives 


which  can  be  solved  for  the 


FB  cos  0B  -  Fs  sin  6  s  =  0 

Fg  sin  e  B  +  Fs  cos  0  g  =  W 

cushion  forces: 

W  cos  e8/cos  (eB  -  0S) 

W 1  sin  9  g/cos  (6  g  -  0  g) 


(88) 


(89) 


Dimensional  considerations  permit  the  coordinates  of  point  P  to  be  written  as 
functions  of  the  thicknesses  tg  and  tg  of  the  compressed  seat  and  back  cushions, 
respectively. 


*Note  that  the  computation  of  Euler  angles  is  required  only  for  the  three- 
dimensional  occupant  model.  For  the  two-dimensional  model,  the  corresponding 
generalized  coordinates  are  obtained  directly  from  input  data. 
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Figure  27.  Forces  acting  on  occupant 
torso  (level  flight). 


Zp  =  Zs  +  (Rx  +  ts)/coses  +  (Xp  -  X$)  tan 9  $ 


Xp  =  Xs  +  (R:  +  tB)/cose  B  -  (Zp  -  Z$)  tano  Q 
which  can  be  solved  for  and  Z^  to  give 

Xp  -  xs  -  fx  sine  jj  +  f2  cose$ 
zp  =  V  ft  cose  B  *  f2  sine  $ 


where 


fl  =  (R1  +  V/cos  (e8  "  0S} 
f„  =  (K,  +  t„)/cos  (e„  -  ec) 


Since  the  force-deformation  characteristics  of  the  cushions  are  known  from 
input  data,  the  compressed  thicknesses  ts  and  tg  can  be  calculated  from  equa¬ 
tion  (89).  These  values,  when  used  in  equation  (91),  give  the  coordinates 
of  point  P  for  the  first  approximation  of  level  flight.  The  equilibrium 
(zero-load)  lengths  of  the  lap  belt  and  shoulder  belt(s)  are  calculated  for 
the  body  in  this  position. 

I 

Next,  the  aircraft  is  rotated  to  the  attitude  specified  by  the  input  condi¬ 
tions  of  pitch,  roll,  and  yaw.  Nose-up  pitch  will  tend  to  load  the  back  cush¬ 
ion,  and  the  analysis  will  be  the  same  as  that  described  above  for  level  flight, 
except  that  the  W'  vector  in  figure  27  will  have  a  component  In  the  XA  direction. 

Nose-down  pitch,  on  the  other  hand,  will  tend  to  load  the  restraint  system. 

An  iterative  procedure  is  used  to  determine  the  correct  position  for  this  case. 
Referring  to  figure  28,  summing  forces  gives  a  set  of  transcendental  equations 


W'  sine^  -  Fs  sines  -  FL  cos  6  L 


0 


W  cos  Q  A  +  F$  cos  6  s  -  FL  sin  6  L  =  0 


(92) 


where  9^  is  the  angle  between  the  floor  and  the  projection  of  the  lap  belt  on 
the  XA  -  ZA  plane.  An  initial  estimate  to  0  L  is  made  using  the  body  position 
calculated  by  equation  (91)  for  the  level  flight  assumption.  The  angle  is 
defined  according  to 

0L  =  sin-1  [(Zp  -  ZL)/  ^  (Xp  -  XL)2  +  (Zp  -  ZL)2]  (93) 

The  forces  in  the  seat  cushion  Fs  and  the  lap  belt  FL  are  determined  using 
this  value  of  e ^  in  equation  (92).  From  the  input  force-deformation  charac¬ 
teristics  for  the  seat  cushion  and  lap  belt,  the  deformations  S<.  and  6^  are 
calculated.  These  deformations  are  used  to  determine  new  values  of  Xp  and  Zp; 
this  procedure  amounts  to  permitting  the  body  to  further  compress  the  seat 
cushion  and  slide  forward  into  the  lap  belt.  Following  through  the  procedure, 
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Figure  28.  Forces  acting  on  occupant  torso 
(nose-down  attitude). 


the  new  length  for  one  side  of  the  lap  belt  is 


L 


L 


(94) 


where  L^e  is  the  equilibrium  length.  The  new  value  of  Xp  is  given  by 

9  o  1/2 

Xp  =  XL  +  C(Ll  -  l/  -  (Yp  -  Yl)-J  COS0L  (95) 


where  is  one-half  the  hip  breadth  and  Yp  is  the  Y-coordinate  of  the  right 
hip  in  the  aircraft  system.  The  new  value  of  Zp  is  computed  for  the  new  cush 
ion  thickness  t^  using  equation  (90),  which  is  repeated  here  for  continuity: 


zp  =  zs  +  (Rj  +  ts)/cose  s  +  (X p  -  Xs)  tan0  $ 


The  new  occupant  position,  determined  by  equations  (9b)  and  (90)  is  used  in 
equation  (93)  to  recalculate  the  lap  belt  angle  0^,  and  the  procedure  is  re¬ 
peated  until  two  consecutive  values  of  Xp  differ  by  less  than  5  percent.  The 
coordinates  of  the  mass  center  of  segment  1  (X^,  Y^,  z^)  are  then  calculated 
from  Xp,  Yp,  and  zp. 


At  this  point  the  generalized  coordinates  for  the  upper  body  have  been  deter¬ 
mined.  The  final  task  will  be  to  determine  the  coordinates  for  the  legs.  Re¬ 
ferring  to  figure  29,  the  angles  Yg  and  ek  can  be  found  from  simple  geometric 
relationships  among  the  dimensions  shown.  The  Euler  angles  H'g,  03.  and  <l>g  are 
obtained  from  Yg  by  using  equations  (8b)  through  (87),  and  the  corresponding 
coordinates  for  segment  10,  by  symmetry.  The  knee  angles  are  given  by 


a  =  "  -  6k  rn  =  9,  11 


to  complete  the  initialization  of  the  generalized  coordinates. 
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Figure  29.  Leg  position, 
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4.3  PROGRAM  SOLUTION  PROCEDURE 


The  first  operation  in  each  solution  step  includes  the  calculation  of  new 
values  for  the  aircraft  acceleration  components  and  their  subsequent  Integra¬ 
tion  to  obtain  aircraft  velocity  and  displacement  components.  Then  the  matrix 
form  of  the  equations  of  motion,  using  equation  (19) 

[A(q)]  {q }  =  {B(q,q)}  +  {P(q)>  +  (R(q,q) }  +  (Q(q\q) } 
are  set  up  for  solution  and  solved,  as  discussed  below. 

4.3.1  Setup  of  Eauatlons  of  Motion 


The  elements  of  [A],  {B } ,  {P},  and  {R}  are  calculated  using  the  current  values 
of  the  generalized  coordinates  and  velocities.  The  elements  of  {Q},  which  is 
the  vector  of  generalized  external  forces,  are  calculated,  as  discussed  In 
section  2.3.  The  external  forces  depend  on  displacements  of  the  aircraft, 
which  determine  the  motion  of  the  seat,  floor,  and  restraint  system  anchor 
points  relative  to  the  body.  From  these  displacements  new  deflections  of  the 
cushions,  floor,  and  restraint  system  are  calculated. 


4.3.2  Solution  of  Equations  of  Motion 


The  system  of  equations  is  solved  for  the  generalized  accelerations  by  first 
combining  the  vectors  on  the  right-hand  side: 

[A]  { q}  3  {  B '}  (97) 


where 


{B*  }  =  {B  }  +  {p}  +  {R}  +  {Q> 


and  solving  for  {q}  using  Crout  decomposition  followed  by  Forsythe-Moler  el  i  mi  - 


nation. 


The  resulting  set  of  N  second-order  differential  equations  have  the  general 


Qj  "  fj  (t»  q^»  cj 2 . q^j»  q^*  ^2’ •  •  • »  qN^ 

qj  (t=°)  =  qjO  *qj  *  qjO  J  “  ]»  2 .  N 


(90) 


where  N,  the  number  of  degrees  of  freedom,  is  12  for  the  two-dimensional  model 
and  29  for  the  three-dimensional  model. 


These  equations  can  be  rewritten  as  2N  first-order  equations  having  the  general 
form 


yj  =  fj  (t»  y^,  y2 .  q2,',,‘  9^) 

qd  -  yj  (99) 

yj  (f°)  ■  9j0  9j  (t-0)  ■  9 j0 

Numerical  integration  of  this  set  of  equations  is  accomplished,  using  the 
Adams-Moulton  predictor-corrector  method  with  a  variable:  step  size.  This 
method  uses  the  difference  equations 


,(P) 


h 


yj^n+l  =  yj,n  +  ^24  (65fi,n  “  59fj,n-l  +  3 /f j  ,n-2  "9fj,n-3) 


(100) 


as  the  predictor  and 


/(.C)  1 
•j  »n+l 


5fj,n-l  +  fd,n-2* 


(101) 


as  the  corrector.  Starting  values  are  provided  by  the  classical  fourth-order 
Runge-Kutta  method.  Input  data  includes  upper  and  lower  error  bounds  for  the 
solution.  Error  bounds  for  each  variable  are  calculated  and  compared  at  each 

step  with  the  difference  between  the  predicted  value  y^  and  the  corrected 

ic  1  J 

value  yl  .  If  this  difference  exceeds  the  upper  bound  for  any  j,  the  step 

size  is  halved.  If  this  difference  is  less  than  the  minimum  error  bound  for 

all  j  and  for  three  successive  steps,  the  step  size  is  doubled.  Halving  the 

step  size  is  accomplished  by  interpolation  of  past  data,  whereas  doubling  is 

effected  by  alternate  selection  of  past  data.  The  solution  can  be  run  with  a 

fixed  step  size  by  makiny  the  upper  and  lower  error  bounds  prohibitively  large 
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and  small,  respectively,  or  by  using  equal  values  for  the  maximum  and  minimum 
step  size  which  are  also  Included  among  input  data, 

4.4  PROGRAM  OUTPUT 


Output  data  consist  of  ten  blocks  of  Information  that  are  selected  for  print¬ 
ing  by  user  input.  The  data  include  time  histories  of  the  variables,  which 
are  simply  sorted  during  the  solution  at  predetermined  print  intervals  as 
fol lows: 


1.  Occupant  segment  positions  in  the  aircraft  coordinate  system  (X,  Y, 
Z,  pitch,  and  rol 1 ), 

2.  Occupant  segment  velocities  in  the  aircraft  coordinate  system  (X,  Y, 
and  Z). 

3.  Occupant  segment  accelerations  in  the  segment-fixed  coordinate  sys¬ 
tems  (x,  y,  z,  and  resultants)* 

4.  Restraint  system  loads. 

5.  Cushion  loads. 

6.  Aircraft  displacement,  velocity,  and  acceleration. 

7.  Injury  criteria.  : 

8.  Details  of  contact  between  the  occupant  and  the  seat  in  front. 

9.  Seat  structure  nodal  forces. 

10.  Seat  structure  element  stresses. 

Printer  plots  are  provided  for  occupant  segment  accelerations,  restraint  sys¬ 
tem  loads,  and  cushion  loads.  The  option  of  two  different  filters  is  also 
provided  for  the  occupant  segment  accelerations  and  cushion  loads. 

Two  of  the  above  blocks  of  output  data  will  be  discussed  in  further  detail. 

4.4.1  Impact  Prediction 

For  prediction  of  impact  between  the  occupant  and  the  seat  in  front,  26  sur¬ 
faces  are  defined  on  the  body.  These  surfaces  'were  illustrated  in  figure  15, 
and  their  dimensions  discussed  in  section  2.5. 
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The  distance  between  each  of  these  occupant  contact  surfaces  and  the  surfaces 
of  the  seat  In  front  is  calculated  during  execution  of  the  program.  When  con¬ 
tact  occurs  between  an  occupant  surface  and  a  contact  surface,  the  time  and 
relative  velocity  of  impact  are  computed  and  stored  for  printing.  The  impact 
conditions  determined  in  this  way  can  be  used  in  evaluation  of  injury  potential 
for  a  given  cockpit  configuration. 

4.4.2  Injury  Criteria 

The  injury  criteria  used  in  the  program  were  selected  as  the  most  suitable  for 
aircraft  crash  analysis.  SOM-TA  output  can  bo  used  to  determine  the  potential 
for  injury  to  four  regions  of  the  body:  head,  vertebral  column,  thorax,  and 
legs.  Each  of  these  is  discussed  in  the  remainder  of  this  chapter. 

4. 4. 2.1  Head  Injury.  An  accepted  criterion  for  head  injury  is  the  Severity 
Index  (SI)  developed  by  C.  W.  Gadd  (references  17  and  18),  which  is  calculated 
for  the  head  and  chest  according  to 


SI  = 


an  dt 


where  a  -  acceleration  in  G  as  a  function  of  time 
n  =  weighting  factor,  2.5  for  head  impacts 
t  =  time  in  seconds. 


(102) 


Although  Gadd  used  uniaxial  acceleration  in  his  validation  of  the  Severity 
Index,  Federal  Motor  Vehicle  Safety  Standard  208  requires  the  use  of  resultant 
acceleration.  A  tolerable  $1  value  of  1000  is  accepted  for  frontal  impact  of 
the  head,  and  a  value  of  1500,  for  distributed,  or  noncontact  accelerations 
(reference  19). 


Also,  the  Head  injury  Criterion  (HIC)  contained  in  Federal  Motor  Vehicle  Safety 
Standard  208  is  calculated  according  to 


(103) 


HIC  =  max  fr — ^—r~  S  adt  J  (t5  -  t.) 

\  2 "  h  ^  1 

where  a  is  the  resultant  head  acceleration  in  G  and  and  tg  are  any  two 
points  in  time  (sec)  during  the  crash  event.  A  tolerable  value  of  1000  is 
accepted  for  the  HIC. 

4.4. 2. 2  Thoracic  Injury.  Although  chest  deflection  Is  commonly  accepted  as 
the  preferred  tolerance  criterion,  it  is  not  presently  predicted  by  SOM-TA. 

FMVSS  208  currently  specifies  as  acceptable  an  acceleration  pulse  which  ".  .  . 
shall  not  exceed  60  G  except  for  Intervals  whose  cumulative  duration  is  not 
more  than  3  msec,"  where  the  acceleration  is  the  resultant  measured  at  the 
center  of  gravity.  This  acceleration  is  printed  by  SOM-TA.  Previously,  FMVSS 
208  had  applied  a  Severity  Index  to  the  chest  acceleration  pulse.  This  index 
was  calculated  in  exactly  the  same  manner  as  the  head  Severity  Index  discussed 
previously,  and  the  limit  of  1000  was  the  same  as  that  for  the  head.  SOM-TA 
also  prints  this  SI  for  the  thorax. 

4. 4. 2. 3  Vertebral  Injury.  For  determination  of  the  potential  of  ^^tebral 
injury,  the  SOM-TA  two-dimensional  occupant  model  has  beam  elements  In  the 
lumbar  spine  and  neck.  Axial  forces  and  bending  moments  are  output  for  these 
segments.  Tolerable  force  levels  for  the  lumbar  spine  have  not  yet  been 
established.  However,  based  on  cadaver  tests,  references  20  and  21  have  pub¬ 
lished  tolerable  bending  moments  for  the  neck  such  as  1700  in. -lb  in  flexion 
(forward  bending)  and  500  in. -lb  ft-lb  in  extension  (rearward  bending),  re¬ 
spectively. 

The  dynamic  response  index  ( DR  I )  is  also  computed  by  SOM-TA  as  a  measure  of 
the  probability  of  spinal  Injury  due  to  a  vertical  acceleration  parallel  to 
the  spine  (reference  22). 

In  this  model,  the  response  of  the  body  to  acceleration  parallel  to  the  spine 
is  modeled  by  a  single  lumped-mass,  damped-spring  system  as  shown  in  figure  30, 
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,‘3. 


*Dynamic  Response  Index  g  -  386  in. /sec' 


Figure  30.  Model  used  for  prediction  of  spinal 
injury  (from  reference  22). 


or,  In  other  words,  the  total  body  mass  that  acts  on  the  vertebral  column  to 
cause  deformation  is  represented  by  the  single  mass.  In  general,  the  motion 
of  the  system  shown  in  figure  29  obeys  the  relationship 


I 


~  +  26uj  +  u,  2<$  =  Z 

dt  n  dt  wn 


(104) 


The  solution,  the  deflection  6,  is  representative  of  the  deformation  of  the 
spine,  and  the  last  term  of  the  left-hand  side  of  equation  (104),  divided  by 
the  gravitational  acceleration,  is  the  OR  I.  The  properties  used  in  the  model 
were  derived  from  tests  involving  human  subjects  and  cadavers.  For  example, 
the  spring  stiffness  k  was  determined  from  tests  of  human  cadaver  vertebral 
segments;  damping  ratios  were  determined  from  measurements  of  mechanical  im¬ 
pedance  of  human  subjects  during  vibration  and  impact.  The  acceleration  input, 
Z,  is  the  component  of  seat  pan  acceleration  parallel  to  the  vertebral  column. 


Equation  (99)  is  solved  lor  S  simultaneously  with  the  occupant  equations  of 
motion,  using  the  constants  r  -  0.224  and  m  --  ‘>2.9  rad/sec.  The  DR  I  is  then 
calculated  at  each  step  by 


DR  I  =  u>n  'Vq 


(105) 


4. 4. 2.4  Leg  Injury 


FMVSS  208  specifies  a  limit  of  2,250  lb  for  the  axial  compressive  load  in  each 
femur.  This  tolerance  level  is  based  primarily  on  the  work  of  Patrick,  Kroell , 
and  Mertz  (references  23  and  24),  who  tested  embalmed  human  cadavers  in  full- 
scale  sled  experiments.  Other  experiments  that  used  a  stationary  seated  cadaver 
struck  by  a  moving  rigid  impactor  produced  higher  fracture  loads,  but  injury 
patterns  were  not  representative  of  field  accidents.  Viano  and  Khalil  (refer¬ 
ence  25)  demonstrated  that  the  poor  correspondence  between  the  latter  tests  and 
actual  accidents  could  be  be  attributed  to  the  pulse  duration  (on  the  order  of 
TO  msec),  significantly  shorter  than  that  which  occurs  in  real-world  accidents. 

If  simulation  of  contact  between  the  occupants  and  the  seat  back  is  requested, 
SOf'l-TA  prints  (and  plots)  the  time  variation  of  the  axial  load  in  each  femur. 

The  time  and  impact  velocity  for  each  such  contact  is  also  printed,  along  with 
identification  of  seat  back  surface. 


5.0  MODEL  VALIDATION 

Dynamic  tests  of  transport  seats  conducted  by  the  FAA  Civil  Aeromedical  Insti¬ 
tute  (CAMI)  in  1983  for  Simula  were  used  in  validation  of  Program  SOM-TA. 

Three  production  seat  designs  and  various  modifications  to  the  production  de¬ 
signs  were  tested.  The  design  that  was  chosen  for  validation  utilized  a  produc¬ 
tion  seat  modified  by  incorporating  energy-absorption  capability  in  the  longi¬ 
tudinal  direction.  This  seat  produced  more  interesting  results  for  validation 
than  others  because  it  underwent  large  deformations  while  maintaining  structural 
integrity  In  the  dynamic  test. 
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5 •  1  Description  of  Seat  Structure 

Figure  31  shows  the  basic  production  seat.  The  seat  pan  assembly  is  built 
around  aluminum  front  and  rear  tubes  which  are  connected  by  four  aluminum 
supports.  The  front  of  each  support  rests  on  the  top  of  the  front  tube  and  is 
attached  by  two  bolts  running  through  the  support  and  the  tube.  The  rear  of 
each  support  has  a  semicircular  cut-out  that  slips  over  the  rear  tube  and  at¬ 
taches  with  one  bolt.  Three  of  the  supports  have  cantilevered  rear  extensions 
which  provide  a  pivot  and  anchor  point  for  the  arm  rests  and  seat  backs.  The 
seat  cushions  rest  on  perforated  aluminum  sheets  which  are  riveted  to  the  front 
tube  In  front,  and  wrapped  around  a  small  aluminum  tube  which  slips  in  between 
the  cantilevered  extensions  in  the  rear. 
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The  ley  assembly  consists  of  front  and  rear  legs  of  formed  stainless  steel 
sheet,  which  are  connected  diagonally  by  a  square  aluminum  tube.  The  front 
legs  cradle  and  bolt  to  the  front  tube,  while  the  rear  legs  art:  held  to  the 
rear  tube  by  steel  straps  which  are  wrapped  around  the  tube  and  bolted  to  the 
legs.  Each  rear  leg  is  connected  to  the  diagonal  tube  by  a  steel  fitting  which 
is  riveted  to  the  leg  and  the  tube,  and  is  also  bolted  to  the  steel  strap. 

The  seat  pan  and  leg  assembly  are  shown  in  figure  32. 

The  rear  legs  are  attached  to  the  floor  track  by  a  fitting  consisting  of  two 
studs  and  a  lever-operated  plunger.  The  plunger  moves  in  and  out  of  the  track 
and  is  in  front  of  the  fitting.  The  fitting  is  attached  to  the  ley  by  a  single 
bolt  and  is  shown  in  figure  33. 
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Figure  33.  Rear  track  fitting  of  a  production 
aircraft  seat. 


The  front  legs  are  attached  to  the  track  by  a  single  stud  which  is  incorporated 
in  an  antirattle  mechanism  shown  in  figure  34.  The  stud  is  bolted  to  the  yoke, 
which  is  attached  to  both  the  front  leg  and  diagonal  tube  by  a  single  bolt.  An 
aluminum  strap  connects  the  front  stud  and  rear  fitting  in  order  to  transfer 
the  loads  from  the  tubular  diagonal  brace  to  the  plunyer. 


The  seat  structure  includes  attachment  fittings  for  standard  lap  belts.  These 
fittings  are  bolted  to  the  forged  aluminum  supports.  !3oth  the  back  and  bottom 
cushions  are  fabricated  of  plastic  foam.  The  back  cushion  is  contoured  for 
comfort  while  the  bottom  cushion  is  flat  on  top  but  curved  on  the  bottom  to  con 
form  to  the  perforated  aluminum  seat  pan.  Each  seat  has  a  food  tray  mounted 
on  its  back,  and  the  seat  back  has  an  adjustable  reclining  mechanims.  The 
two  middle  armrests  can  pivot  up  flush  with  the  seat  backs. 


Figure  34.  Front  track  fitting  of  a  production 
aircraft  seat. 

The  production  seat  design  was  modified  at  Simula  to  provide  longitudinal 
energy-absorbing  capability  by  replacing  the  diagonal  brace  with  a  compressive 
Inversion  tube  energy  absorber.  A  sketch  of  the  seat  kinematics  using  the 
compressive  energy  absorber  Is  shown  In  figure  35.  In  order  for  the  seat  to 
undergo  the  approximate  curvilinear  translation  the  sketch  depicts,  the  front 
and  rear  legs  needed  to  be  released  to  act  as  a  linkage. 

The  energy  absorbers  used  are  shown  in  figure  36.  Both  employ  the  load  atten¬ 
uation  method  of  inverting  an  aluminum  tube.  Due  to  the  asymmetry  of  the  seat 
the  window-side  energy  absorber  is  required  to  stroke  at  twice  the  load  of 
the  aisle-side  energy  absorber,  and  thus  uses  two  inversion  tubes.  The  window 
side  and  aisle-side  energy  absorbers  begin  to  stroke  at  2800  lb  and  1400  lb, 
respectively. 

The  track  fittings  were  also  modified  to  provide  release  about  pitch  and  roll 
axes.  Both  front  and  rear  track  fittings  had  three  studs  and  were  locked  in 
both  positions.  The  modified  track  fitting  is  shown  in  figure  37. 


t-'igure  37.  Modified  transport  aircraft 
track  fittings. 


mamic  Test  Conditions 


The  seat  design  was  tested  at  CAMI  in  the  forward-facing  configuration  with  a 
9-G,  50  ft/sec.  sled  deceleration  (figure  38).  Three  50th-percentile  anthro¬ 
pomorphic  dummies  were  used  to  simulate  the  seat  occupants  (figure  39).  Six- 
axis  load  cells  were  placed  under  the  seat  legs  to  measure  the  seat  reactions 
exerted  on  the  floor.  Lap  Pelt  forces  were  also  measured  on  both  sides  of  the 
lap  belt  for  each  occupant.  Accelerometers  were  placed  at  dummy  pelvis,  chest, 
and  head  locations. 


Prior  to  dynamic  testing,  the  test,  fixture  was  adjusted  to  simulate  floor 
warping  conditions.  The  aisle  side  of  the  seat  was  pitched  down  10  degrees 
(figure  40)  and  the  window  side  of  the  seat  was  rolled  outwards  10  degrees 
( figure  41 ) . 

5.3  SOM-TA  Simulation  Results 

The  finite  element  model  of  the  seat  structure  illustrated  in  figure  4?!  con¬ 
sists  of  26  nodes  and  33  beam  elements.  The  complete  input  data  listing  is 


presented  in  figure  43. 


The  compressive  energy  absorbers  were  modeled  with  beams  21  and  24  between 
nodes  9-12  and  18-21,  respectively.  The  ends  of  the  energy  absorbers  were 
modeled  as  pin-socket  joints  by  using  the  appropriate  bean  end  release  codes. 
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Figure  43  (contd) 


Input  data  listing  for  simulation 
of  CAMI  Test  A83-085. 
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Figure  43  (contd).  Input  data  listing  for  simulation 

of  CAMI  Test  A83-085. 


The  nonlinear  force-deflection  characteristics  of  the  energy  absorbers  were 
modeled  using  a  fictitious  material  property  data  for  the  corresponding  axial 
stress  and  tangent  modulus.  All  other  beau  elements  made  of  4130  steel  and 
2024-T3  aluminum  were  modeled  as  elastic  beams. 

The  floor  warping  condition  of  10  degrees  pitch  on  the  aisle  side  was  simu¬ 
lated  by  enforcing  a  downward  displacement  of  2.68  In.  at  node  9.  Track  fit¬ 
ting  design  already  provided  release  for  floor  warping  about  the  roll  axis. 

The  three  (>Oth-percentile  dummies  were  modeled  using  standard  (two-dimensional) 
occupants,  with  occupant  segment  data  parameters  available  in  the  program. 

The  data  set  in  figure  43  requests  all  occupant  output  data  for  the  center  oc¬ 
cupant. 

As  evidenced  by  the  post-test  view  In  figure  44,  the  seat  did  not  experience 
any  failures  in  the  dynamic  test.  The  floor  deformation  prior  to  the  test  did 
not  affect  the  performance  of  the  seat,  due  to  the  releases  provided  by  the 
track  fittings. 

The  predicted  occupant  and  seat  positions  from  SOM-TA  simulations  are  shown  in 
figures  46  and  46,  respectively.  The  occupant  chest  contacts  the  knees  at  about 
200  msec.  Program  SOM-TA  does  model  the  interaction  between  the  chest  arid  the 
knees  using  an  exponential  force-deflection  relationship  and  a  viscous  damper. 

Program  SOM-TA  predictions  for  the  center  occupant  head,  chest,  and  pelvis  ac¬ 
celerations,  lap  belt  forces,  and  left-seat  leg  reactions  are  compared  with 
the  measured  test  data  in  figure  47  through  figure  S3.  Occupant  accelerations 
have  peaks  around  200  msec,  an  indication  of  the  contact  between  the  occupants' 
ciiest  and  knees. 
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PROGRAM  SOM-TA 

TRANSPORT  AIRCRAFT  SEAT  (MODIFIED) 
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TRANSPORT  AIRCRAFT  SEAT  (MODIFIED) 
TIME  -  0.0400  SEC. 
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TRANSPORT  AIRCRAFT  SEAT  (MODIFIED) 
TIME  -  0.3000  SEC. 


PROGRAM  SOM-TA 

TRANSPORT  AIRCRAFT  SEAT  (MODIFIED) 
TIME  =  0.3000  SEC. 


Figure  46.  Program  SOM-TA-predicted  seat  displacements  for  modified 
transport  aircraft  seat  ( CAM  I  Test  A83-085). 
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6.0  CONCLUSIONS 


The  Seat/Occupant  Model  -  Transport  Aircraft  (Program  SOM-TA)  provides  a  method 
for  analyzing  the  performance  of  a  seat  and  occupant  in  an  aircraft  crash  en¬ 
vironment.  For  preliminary  analyses  where  details  of  the  seat  structure  are  not 
known,  a  rigid  seat  can  be  used.  Otherwise,  finite  element  modeling  capability 
is  available  to  model  the  seat  structure  in  greater  detail. 
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APPENDIX  A 

OCCUPANT  SEGMENT  POSITION: 
THREE-DIMENSIONAL  MODEL 


Referring  to  figure  A-l,  and  using  pn  *  Ln  -  pn  the  absolute  position  of  the 
mass  center  of  each  body  segment  Is  given  below.  The  elements  of  the  trans¬ 
formation  matrices  [Tn]  are  functions  of  the  generalized  coordinates,  as  given 
by  equations  (6),  (9),  and  (10). 

Segment  1: 


(Xj,  Y^,  Zj),  the  coordinates  of  the  reference  point  on  the  body  are  the  gen¬ 
eralized  coordinates  (q ^ ,  q2,  q3). 


Segment  2: 


Segment  3: 


Segment  4: 


igure  A-i.  Body  segment  lengths  and  mass  center 
locations  for  three-dimensional  model. 
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APPENDIX  B 


OCCUPANT  SEGMENT  POSITION 
TWO-DIMENSIONAL  MODEL 


Referring  to  figure  6,  the  absolute  position  of  the  mass  center  of  each  body 
segment  is  given  below.  L$  and  are  the  lateral  distances  from  the  mid- 
saggital  plane  to  the  shoulder  and  hip  joints,  respectively. 

Segmt'  it  1 : 

(Xj,  Yp  Zl),  the  coordinates  of  the  reference  point  on  the  body.  Xj^  and  Z ^ 
are  the  generalized  coordinates  q1  and  q^. 
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Director  DuCentre  Exp  DE  LA 

(1) 

Trans portat ion 

i  Center  Library 

Navigation  Aerineene 

Evanston,  ILL 

60201 

9141  Orly,  France 

ANE-40 

(2) 

ACT- 62 4 

(2) 

ASW-53B 

(2) 

ASW-52C4 

(2) 

AAL-62 

(2) 

AAC-44.4 

(2) 

APM-13  Nigro 

(2) 

M-493.2 

(5) 

ACE-66 

(2) 

Bldg.  10A 

AEA-66.1 

(3) 

ADL-1 

(1) 

APM-l 

(1) 

ADL-32  North 

(1) 

ALG-30O 

(1) 

APA-300 

(1) 

AES-3 

(1) 

ACT- 8 

(1) 

AGL-60 

(2) 

ANM-60 

(2) 

FAA,  Chief,  Civil  Aviation  Assistance  Group  (1) 

Al  Astorga 

(1) 

Madrid,  Spain 

Federal  Aviation 

c/o  American  Embassy 

Administration  (CAAG) 

APO-Naw  York  09285-0001 

American  Embassy,  Box 

38 

APO-  New  York  09285-0001 


Dick  Tobiason  (1) 

ATA  of  America 

1709  New  York  Avenue,  NW 

Washington,  DC  20006 


Government  Activities 

FAA,  Washington,  DC  20591 
(Attn:  Harold  W.  Becker,  ASF-300; 

Thomas  McSweeny,  AWS-100)  (2) 

FAA,  4344  Donald  Douglas  Drive, 

Long  Beach,  CA  90808 

(Attn:  Stephen  Soltis,  ANW-102N)  (1) 

FAA,  Central  Region  Headquarters, 

601  East  12th  Street,  Federal  Building, 

Kansas  City,  Missouri,  64106 

(Attn:  Earsa  Tunkesley,  ACE-100)  (1) 

FAA,  Southwest  Region  Headquarters, 

4400  Blue  Mound  Road, 

P.0.  Box  1689 

Fort  Worth,  Texas  76101 

(Attn::  John  Shapley,  ASW-110)  (1) 

FAA,  Northwest  Mountain  Region 
Head  quarters , 

17900  Pacific  Highway  South,  C-68966, 

Seattle,  Washington  98168 

(Attn:  Steven  Wallace,  ANW-110)  (1) 

FAA,  Mike  Monroney  Aeronautical  Center, 
P.O.  Box  26082, 

Oaklonoina  City,  OK  73125 

(Attn:  Ricn.ird  Chandler,  AM-119)  (1) 

NASA,  Langley  Research  Center, 

Hampton,  VA  2 336 5 

(Attn:  Emilio  Alfaro~Bou,  MA-493; 

Huey  Cardin,  MS-495)  (2) 

U.S.  Army  Aviation  Applies  Technology 
Directorate,  USAAKTA,  (AVSCOM) , 

Fort  Eustis,  VA  2360  4 

(Attn:  Roy  Burrows,  Code  SAVRT-TY-ASV) 

U.S.  Navy  Naval  Air  Development  Center, 
Warminister,  PA  18974 

(Attn:  Leon  Domzalski,  Code  60322)  (1) 

NTSB,  800  Independence  Ave .  S.E., 

Washington,  DC  20594 

(Attn:  John  C.  Clark,  TE  60)  (1) 


Non-Government  Ac  t  iv it ies 

Beech  Aircraft  Corp., 

P.O.  Box  85, 

Wichita,  KS  67201 
(Attn:  Dayton  L.  Hartley; 

James  E.  Terry;  William  Schultz)  (3) 

Bell  Helicopter  Co., 

P.O.  Box  482, 

Fort  Worth,  TX  76101 

(Attn:  James  Cronkite,  Dept,  81, 

MS  11;  Roy  G.  Fox)  (2) 

Boeing  Airplane  Co., 

P.O.  Box  3707, 

Seattle,  WA  98124 

(Attn:  Edward  Widmayer,  MC-9W-22)  (1) 

Boeing  Co.,  Vertol  Division, 

P.O.  Box  16858, 

Philadelphia,  PA  91942 

(Attn:  Denise  Vassilakos,  MS  P30-27)  (1) 

Cessna  Aircraft  Co., 

P.O.  Box  7704, 

Wich  ita ,  KS  672  77 

(Attn:  John  Berwick;  Robert  Held; 

Richard  Soloski)  (3) 

Fairchild  Aircraft  Corp., 

P  .0  .  Box  3246 , 

San  Antonio,  TX  78284 
(Attn:  Walt  Dwyer)  (1) 

General  Dynamics /Conva ir , 

P.O.  Box  80847, 

San  Diego,  CA  92138 

(Attn:  L.  Mastay,  MA  80-6030)  (1) 

(2)  Grumman  Aerospace  Corp., 

So.  Oyster  Bay  Road,,  Bethpage, 

L .  I .  ,  NY  11714 
(Attn:  Robert  Winter,  A08-35; 

Allan  B.  Oifko,  A08-35)  (2) 

Gulfstream  Aerospace  Corp., 

P.O.  Box  2206, 

Savannah,  GA  31402 
(Attn:  George  Westphal)  (1) 
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Non-Government  Activities  (Continued) 


Gulfstream  Aerospace  Corp., 
P.0.  Box  22500, 

Oklhoma  City,  OK  73123 
(Attn:  Richard  Southard)  (1) 

Lockheed-Cal iforn ia  Co., 

Burbank,  CA  91503 

(Attn:  Gil  Wittlin,  D  76-12, 

B  63G ,  PLT  A-l)  (1) 

McDonnell  Douglas  Corp., 

3655  Lakewood  Drive, 

Long  Beach,  CA  90646 
(Attn:  J.  Webster; 

John  L.  Galligher)  (2) 


McDonnell  Douglas  Helicopter, 
4645  S .  Ash  Ave . , 

Tempe ,  AZ  65182 
(Attn:  Lyndon  Landborne; 

J.K.  Sen)  (2) 

Piper  Aircraft  Corp., 

2925  Piper  Drive, 

Vero  Beach,  FL  32960 
(Attn:  Marion  Dees)  (J. ) 

Sikorsky  Aircraft, 

North  Maint  Street, 

Stratford,  CT  06601 

(Attn:  Brain  Cornell,  MS  5207A; 

Pramonik  Mukunda,  MS  5207A)  (2) 
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